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Summary 

This paper describes the design, analysis, and 
nonlinear simulation results (batch and piloted) for 
a longitudinal controller scheduled to be flight-tested 
on the High- Alpha Research Vehicle (HARV). The 
HARV is an F-18 airplane modified for and equipped 
with multiaxis thrust vectoring. The longitudinal 
controller consists of feedback and feed-forward con- 
trollers. The paper includes a description of the fa- 
cilities, a detailed review of the feedback controller 
design, linear analysis results of the feedback con- 
troller, a description of the feed- forward controller 
design, nonlinear batch simulation results, and pi- 
loted simulation results. 

Linear analysis shows that the feedback controller 
meets the stability design guidelines, is highly robust 
to stability and control derivative parameter changes, 
has good closed-loop damping ratios, and attenuates 
the structural modes well below the design guideline. 

Batch simulation results include maximum pitch 
stick agility responses, angle-of-attack a captures, 
and a regulation for full lateral stick rolls at several 
a’s. The controller meets most of the agility guide- 
lines throughout the a range, has good a regulation, 
and has consistently good transient response for all 
altitudes within the design region. 

Piloted simulation results include task descrip- 
tions for several types of maneuvers, task guidelines, 
the corresponding Cooper-Harper ratings from three 
test pilots, and some pilot comments. The ratings 
show that desirable criteria were achieved for almost 
all of the piloted simulation tasks, most of which 
had very stringent criteria. The control system was 
shown to be extremely departure resistant. 

Introduction 

Development of the short-to-medium range aN- 
aspect missiles have increased the need for a fighter 
airplane to rapidly point its nose at an opponent 
to obtain first shot opportunity during air combat 
maneuvering. Rapid nose pointing is of particular 
importance within visual range. Even though long- 
range missiles are also continually improving, fail- 
ures and rules of engagement will inevitably produce 
close-in encounters. During close-in air combat, the 
need to acquire a position in the opponent’s rear 
quarter has diminished with the continual improve- 
ment of short-to-medium range all-aspect missiles. 
Due to the enhanced capability of these missiles, a 
great deal of research has gone into developing high- 
performance airplanes that can rapidly maneuver at 
angles of attack a greater than the maximum lift, 


which generates very high nose-turning rates. Simu- 
lation studies have shown that operational capability 
and good controllability in this high-a flight regime 
provide tactical maneuver advantages during close- 
in air combat (refs. 1-4) with current missiles and 
guns. The success in simulated air combat studies 
have lead to research airplanes that have the capa- 
bility to maneuver at poststall angles of attack. 

In recent years, four high- a flight test airplanes 
(i.e., X-29, X-31, the High- Alpha Research Vehicle 
(HARV), and the F-16 Multi- Axis Thrust Vector- 
ing (MATV)) have been maneuvering successfully at 
poststall angles of attack which was not previously 
achievable. (See ref. 5.) The X-29 airplane can main- 
tain controlled flight up to a = 40°; the X-31 airplane 
and HARV are designed to maneuver at a’s up to 70° 
with the application of multiaxis thrust vectoring. 
The F-16 MATV airplane has done sustained ma- 
neuvers up to a = 84° with excursions beyond that. 
Thrust- vectored controls provide the X-31, HARV, 
and F-16 MATV with the added control power re- 
quired to maneuver in the high-a flight regime which 
is not achievable with purely aerodynamic controls. 

The HARV is a modified F-18 airplane that began 
flight evaluation in midyear 1991. The main modifi- 
cation is the addition of thrust- vectored controls in 
the pitch and yaw axes, which resulted in an airplane 
that is approximately 4000 lb heavier than an un- 
modified F-18 airplane. The original thrust-vectored 
control laws were developed jointly by NASA and 
McDonnell Aircraft Company and are being used 
during the initial HARV flight evaluations that in- 
clude flight envelope expansion and maneuvering ca- 
pability. Revised controllers based upon modern con- 
trol design methodologies and processes (ref. 6) will 
then be installed on the HARV after the initial flights 
are completed. One of these revised controllers is the 
longitudinal controller described in this paper. 

For completeness, much of the information in ref- 
erence 7 which pertains to the longitudinal controller 
design is also contained in this paper. In particular, 
reference 7 contains the following: 

1. An overview of the feedback controller with 
the application of variable-gain output feedback 
which generates a feedback gain functional that 
incorporates optimal gain scheduling 

2. A derivation and description of the feed-forward 
command generator (FFCG) for pilot stick sensi- 
tivity and mode selection 

3. Formulation of the command generator tracker 
(CGT), which is a feed- forward control 

4. Control system design guidelines 



5. Some linear analysis and nonlinear batch simula- 
tion analysis 

A key part of reference 7 is the design procedure 
for selection and adjustment of the optimal weights 
used in the variable-gain output feedback design al- 
gorithm. Note that the variable-gain methodology 
(refs. 7-10) used in the feedback design is an inte- 
grated approach for the simultaneous design of the 
many operating conditions within the flight regime. 

This paper contains a detailed description of the 
longitudinal controller design, analyses, and non- 
linear simulation (batch and piloted) results. The 
longitudinal controller is comprised of a feedback 
controller and a feed-forward controller. The feed- 
back controller design section includes a description 
of the modeling and optimal weighting formulation 
as well as tables of all design weights. The feedback 
controller is the same design as that described in ref- 
erence 7 but more specific information is presented 
in this paper. 

Linear analysis results for the feedback controller 
design are described in another section of this pa- 
per. Results for 39 design conditions are shown when 
practical; however, three example cases (i.e., low -a 
high-Mach, medium-a, and high-a) are selected to 
show results for those analyses that require an exten- 
sive amount of data or plots. Linear analysis results 
include closed-loop eigenvalues, stability margins, ro- 
bustness to changes in various parameters, and the 
loop transfer function for a combination rigid-body 
and servoelastic frequency response. Tables are in- 
cluded to show all variable-gain feedback components 
and final gains for the three example cases. 

The FFCG configuration described in reference 7 
has been modified. Key equations for the two com- 
mand modes of a and load factor have not changed 
but the implementation is revised to eliminate prob- 
lems with mode transition. The main implementa- 
tion change relates to the a-trim bias. The modifica- 
tions, the key equations as derived in reference 7, and 
implementation are described. A major change is in 
the CGT, which was described as part of the feed- 
back controller in reference 7. The CGT gain was 
too sensitive at several flight conditions and variable 
limits had to be incorporated. The new configuration 
is described and shown in this paper. 

Results for both six-degrees-of-freedom non- 
linear batch and piloted simulations are included 
in this paper. Batch simulation includes selected 
agility responses for comparison with design guide- 
lines (ref. 11) as well as several capture and reg- 
ulation maneuvers to illustrate performance capa- 
bility. The piloted simulation section includes task 


scenarios, task guidelines, and Cooper-Harper rat- 
ings (CHR) for the various maneuvers. Piloted tasks 
were developed to evaluate the longitudinal control 
law acquisition capability, tracking capability, and 
resistance to coupling (regulation of a during rolling 
maneuvers). Cooper-Harper ratings and pilot com- 
ments are given for the version of the control law to 
be flight tested. Some of the problems encountered 
with earlier versions during piloted simulations are 
discussed. 

Symbols 

In the following list of symbols, matrices and 
vectors are shown in boldface; scalars are shown in 
italics. 

A continuous state matrix 

A m discrete weighting matrix for cross 

product of states and controls 

A p continuous plant state matrix 

B continuous control matrix 

B p continuous plant control matrix 

B w continuous plant process noise 

disturbance matrix 

B continuous process noise distur- 

bance matrix 

C plant and controller state to output 

matrix 

Cj steady-state- normalized coefficients, 

i = 1,2,. or 3 

C p plant state to output matrix 

Ci, C 2 ,C 3 steady-state- normalized coefficients 

for a, q , and n Zl respectively 

D plant control to output matrix 

E expectation operator 

El command- generator-tracker feed- 

forward gain 

e exponential function 

/ scalar weighting for cost function 

g gravity, ft /sec 2 

B zx matrix relating plant states to 

integrator states 

H zy matrix relating measurements to 

integrator states 

I W ,I ^ identity matrices of appropriate 

dimension 
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Ixxi lyyi Izz 

moments of inertia about 

Q 

the X-, y-, and Z-axes, 
respectively, slug-ft 2 

R 

Ry 

I'XZ 

product of inertia in 


x z frame, slug-ft 2 

r 

J 

global cost 

s 

Jc 

continuous local cost 

L 

7/^ 

Jd 

discrete local cost 

Urp 

V 

K 

feedback-gain matrix 

V 

K nz 

proportional feedback gain for 
n z , (deg/sec )/g 

^0 

K q 

proportional feedback gain for q 

V C 

Ku 

control filter feedback gain, sec -1 

V 

K „ 

proportional feedback-gain matrix 

^0 

K z 

integrator feedback gain 

w 

K a 

proportional feedback gain for 
a, sec" 1 

Wo 

K 0 ,K 4. 

variable-gain feedback matrix 
partitions 

w 

k 

sample time integer 

w p 

M 

number of operating conditions 

wo 

N 

integer for series summation 

^acc 

n z 

load factor, g units 

X 0 

n Z,C 

load factor command, g units 

p 

body-axis roll rate 

X 

Ps 

static pressure, lb/ft 2 

X 

p 

gain-schedule parameter 

Xa 

Q 

discrete state-weighting matrix 

X C g 

Q 

continuous output-weighting matrix 

x m 

Qc 

impact pressure, lb/ft 2 

Xp 

Qn z 

continuous weight for n z 


Qq 

continuous weight for q 

X q 

Qu 

continuous weight for control filter 

x s 

Q y 

continuous-weighting matrix for 

X u 

outputs 

X u 

Qz 

continuous weight for integrator 

X v 

Qa 

continuous weight for a 

x w 

Q 

pitch rate, deg/sec 

X Z 


pitch acceleration, deg/sec 2 
discrete weight for control 
continuous weight for control 
body-axis yaw rate 
Laplace transform variable 
time, sec 

plant control input vector 

total airspeed, ft /sec 

covariance matrix for sensor noise 

covariance for control input 
pseudonoise 

rate command, deg/sec 

sensor noise vector 

control input pseudonoise 

total process noise covariance 
matrix 

covariance matrix for external 
disturbance 

process noise vector 

plant process noise vector 

external disturbance noise vector 

distance along X-axis from c.g. to 
accelerometer 

state covariance matrix for initial 
condition error 

plant and controller state vector 

time derivative of x 

actuator state vector 

c.g. location along the X-axis, 
percent of mean aerodynamic chord 

longitudinal state vector 

plant state vector 

first derivative of plant state vector 

pitch rate state 

sensor state vector 

control position state 

derivative of control position state 

airspeed state 

washout filter state 

integrator state 
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X Z 

Xa 

X0 

x 0 

y 

Vcmd 

Vn z 

y p 

Vq 

Vu 

Vz 

Va 

z 

z cg 

a 

a 


&OC 

a oc,lg 


(3 

r 

T p 

Tty 


'Jw 

An ZjC 

AT 

Aoc oc 

6s 

6 V 

6ij 

6sc 


time derivative of x z 

angle-of-attack state 

pitch attitude state 

state vector for initial condition 
errors 

plant and controller output vector 

command from feed-forward com- 
mand generator 

output for load factor 

plant output vector 

output for pitch rate 

output for controller position 
command state 

output for integrator state 

output for angle of attack 

^-transform variable 

c.g. location along Z-axis, relative 
to waterline (WL) 

angle of attack, deg 

time derivative angle of 
attack, deg/sec 

angle-of-attack command, deg 

angle-of-attack command trim, deg 

angle-of-attack command trim for 

1 9 , deg 

sideslip angle, rad 

discrete plant and controller control 
matrix 

discrete plant control matrix 

discrete plant and controller process 
noise matrix 

discrete plant process noise 
perturbation in n^ c , g units 
sampling period, 0.0125 sec 
incremental bias trim command 
metric stabilator control, deg 
pitch thrust vector control, deg 
Kronecker delta 
stabilator command, deg 


fisp 

pilot stick command, in. 

five 

pitch thrust- vectoring 
command, deg 

V 

plant and controller measurement 
noise vector 

Vp 

plant measurement noise vector 

e 

pitch attitude, deg 

M 

structured singular value 

V 

measured variables used to calcu- 
late p 


discrete plant and controller state 
transition matrix 

4> 

bank angle, deg 


discrete plant state transition 
matrix 

Subscripts: 


c 

controller or command 

i,j 

series integers 

k 

coefficient for sampling sequence 

P 

plant 

Superscripts: 


T 

transpose 

-1 

inverse 

Special notation: 

*(•••) 

arguments of variable X 

Abbreviations: 


A/B 

afterburner (maximum thrust) 

CGI 

computer-generated image 

CGT 

command generator tracker 

CHR 

Cooper-Harper rating 

c.g. 

center of gravity 

DMS 

Differential Maneuvering Simulator 

FFCG 

feed-forward command generator 

HARV 

High- Alpha Research Vehicle 

HUD 

heads-up display 

KCAS 

knots, calibrated airspeed 

PIF 

proportional integral filter 

PIO 

pilot-induced oscillation 

RFCS 

Research Flight Control System 
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Description of Facilities 

Airplane 

The configuration used for this control law de- 
sign is an F-18 airplane modified to have multiaxis 
thrust vectoring for additional pitch and yaw con- 
trol power. The basic F-18 airplane is a multirole 
fighter-attack airplane with supersonic dash capa- 
bility and, by current standards, good low-speed, 
high-a: maneuverability. Thrust- vectoring capability 
was achieved by removing the secondary (divergent) 
nozzles and adding three thrust- vectoring vanes per 
engine. Major dimensions and key features of the 
F-18 airplane with thrust vectoring are shown in fig- 
ure 1. The modified airplane, known as the High- 
Alpha Research Vehicle (HARV), is discussed in ref- 
erence 12, and is shown in figure 2. 

The simulation normally used HARV weight and 
inertias corresponding to 60-percent internal fuel. 
Heavy (maximum fuel) and light (minimum landing 
fuel) configurations were also tested; however, only 
results from the nominal (60-percent internal fuel) 
configuration are given in this report. Weights, 
inertias, and center-of-gravity (c.g.) locations for all 
configurations are shown in table I. 

Simulation 

The HARV simulation was built from the pro- 
duction F-18 airplane nonlinear aerodynamic, engine, 
and control system models obtained from McDonnell 
Aircraft Company (McAir). The F-18 airplane simu- 
lation on which the HARV model is based is discussed 
in detail in reference 13. The HARV thrust-vectoring 
capabilities are discussed in reference 14. The origi- 
nal McAir aerodynamic database covers the a range 
from —10° to 90°, the sideslip angle /3 range from 
—20° to 20°, altitudes up to 60 000 ft, and speeds up 
to Mach 2.0. The F-18 airplane simulation was mod- 
ified to account for the thrust- vectored engine and its 
effect on the configuration aerodynamics. 

The HARV engine model (obtained from McAir 
who also designed the thrust- vectoring hardware for 
the airplane) incorporated thrust-vectoring capabil- 
ity, accounted for the effects of Mach number and 
altitude, and included the dynamic response of en- 
gine thrust. The effects of a and vane deflection on 
thrust were also included. Gross thrust and ram drag 
were tabulated separately allowing thrust vectoring 
to affect gross thrust only. The effect of a on engine 
thrust was obtained with a General Electric engine 
simulation having a high level of fidelity and the in- 
let rake survey data from F-18 airplane flight tests. 


Effects of vane deflection on gross thrust and jet de- 
flection were obtained from ground-based tests. (See 
ref. 14.) 

Aerodynamic increments were added to the data- 
base to account for the addition of thrust- vectoring 
vanes, actuator housings, and spin chute. Jet- 
induced effects were added to account for the change 
in airflow over the airframe due to thrust vectoring. 
This data was obtained from wind tunnel tests de- 
scribed in reference 15. 

A Research Flight Control System (RFCS), which 
consisted of a new lateral-directional control law 
and the longitudinal control law discussed herein, 
replaced the existing F-18 airplane control system. 
The thrust- vectored commands from the RFCS go to 
a vane control system known as the mixer-predictor. 
The mixer-predictor converts pitch, yaw, and roll 
thrust-vectoring commands into commands for the 
six thrust- vectoring vanes to yield the required jet 
deflection. 

Piloted Simulator 

The piloted simulation portion of the control law 
design process was conducted with the math model 
described above in the NASA Langley Differential 
Maneuvering Simulator (DMS). The DMS is a fixed- 
base simulator that has the capability of simulat- 
ing two airplanes as they maneuver relative to each 
other and the Earth. A wide-angle visual display is 
provided for each pilot. As illustrated in figure 3, 
the DMS consists of two 40-ft-diameter projection 
spheres each enclosing a cockpit, airplane image pro- 
jection system, and computer-generated image (CGI) 
sky-Earth-Sun projection system. Each pilot is pro- 
vided with a projected image of the opposing air- 
plane, which gives range, position, and orientation 
cues. Reference 16 contains a detailed, although not 
current, description of the DMS. 

The DMS is driven by a real-time digital sim- 
ulation system built around a CONVEX 3800-series 
computer 1 . The dynamics of the airplane and control 
system were calculated using six- degrees-of- freedom 
rigid-body equations of motion with an 80-Hz frame 
rate. Input and output to and from the computer, 
the cockpit, and the displays had a 40-Hz frame 
rate. Overall transport delay of the system is around 
110 msec. 

A photograph of a cockpit and target visual dis- 
play is shown in figure 4. Each cockpit incorporates 


Convex Computer Corp., Richardson, TX. 
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three CRT heads-down displays and a heads-up dis- 
play (HUD) with a computer-driven gunsight repre- 
sentative of current fighter aircraft equipment. For 
this study, a fixed reticle projected on the heads-up 
display was used for tracking. The simulator displays 
are similar to F-18 airplane displays with some mi- 
nor modifications to facilitate some of the maneuvers 
and the tracking tasks. 

A movable center stick was provided for pitch and 
roll commands from the pilot. Longitudinal stick 
forces and gradients were configured to model those 
of the F-18 airplane. Longitudinal stick travel was 
2.5 in. forward and 5.0 in. aft with a force gradient 
of 7 lb/in. and a breakout force of 2 lb. 

Feedback Controller Design 

This section contains the mathematical equations 
of the HARV feedback control design. The argu- 
ment p shown in the equations represents one or 
more of the parameters that change during flight. 
Optimal gain-scheduling is accomplished through the 
argument p. 

For all equations within this report, matrices 
and vectors are shown in boldface. In the specific 
design example where a matrix or vector is used as a 
scalar, the symbol will be shown in italics. Equations 
presented in reference 7 are general; however, the 
equations presented in this paper are related to the 
specific controller and, therefore, some of the symbols 
will be shown as scalars. 

Modeling 


where v c is the rate command (control feedback 
point), x u is the control position state, x z is the 
integrator state, and H zy and H zx are matrices 
that select the feedback measurements and their 
equivalent states to be integrated. Because there is 
a single control, only one integrator state is used. 

The combination of plant and controller dynamics 
gives the following general equation 

x(p) = A(p)x(p) + B(p)m(p) + B(p)w(p) 1 

? ( 4 ) 

y(p) = c(p)x(p) + 77(p) J 


which can be expanded as 


r x p (p) 

< *«(p) 
l ®*(p) 


" A p (p) 

B P (p) 

O' 

( x p(p) ) 

01 

0 

0 

0 

< *u(p) > + 

In 

,H**(p) 

0 

0. 

l s*(p) J 

. 0 _ 


Vc(p) 



(5a) 


y P (p) 

Vu(p) 
Vz{ p) 


'Cp(p) 0 O' 
0 I u 0 
0 0 1 *. 


Xp(p) 

*u(p) 

Mp) 



(5b) 


in which the terms I u and I represent scalars with 
magnitude one. In general I u and I z are identity 
matrices of appropriate dimension. 

The plant state x p is partitioned into four states 
as 


Continuous model. The dynamic process for 
the plant is represented by 

xp(p) = Ap(p)x p (p) + Bp(p)w p (p) + B ro (p)wp(p) 

y P (p) = c P (p) x P (p) + »fr(p) 

in which x ; , and y p are the plant state and output 
vectors; u p is the scalar plant control; w p and rj p 
are plant process and measurement noise vectors; 
A p , Bp, and C p are the plant state, control, and 
output matrices, respectively; and B w is the plant 
process noise disturbance matrix. Each noise process 
is assumed to be white with zero mean; the processes 
are uncorrelated with each other. 

Controller dynamics are also described by state 
equations 

& u (p) = «c(p) ( 2 ) 

£*( p) = H* w (p)y p (p) = H 2 X (p)x p (p) + H 2J/ (p)? 7 p(p) (3) 



Xp( P)=[x£(p) x£(p) xj’(p) x w (p)] T (6) 

For the HARV longitudinal control law, x m repre- 
sents four longitudinal states, x a represents six actu- 
ator states, x s represents four sensor (or filter) states, 
and x w is a washout filter state. The washout filter 
is actually part of the controller. However, for design 
and some types of analysis, the filter is included at 
the input to the plant because the controller output 
splits into two signals; one is a direct output signal 
for the stabilator command and the other is an out- 
put signal to the washout filter. The combination of 
the washout filter and the plant allows design with a 
single control (dimensions of x u and v c ). The total 
dimension of x p is 15 states and the total dimension 
of x is 17 states because of the inclusion of one con- 
trol and one integrator. The four longitudinal states 
used in the HARV design are 

Xm(p) =[^(P) ^a(p) x q(p) ^(P)] T ( 7 ) 
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which represent total airspeed state x v , angle-of- 
attack state x a , pitch rate state x q , and pitch at- 
titude state xq. 

The plant output y p has three components, which 
represent measurements of angle of attack a , pitch 
rate q , and load factor n*, and is shown as 

yp(p)=[Va(p) Vq{ P) J/n 2 (p)] T (8) 

which means that the dimension of C p is 3 x 15. 


Discrete model. The design approach is to 
discretize the continuous equations (plant plus con- 
troller) to form 


{ x P (p, k + 1) 
x u (p,k + i) 
x z (p,k + 1) 


*j>(p) 

0 

,(AT)H«(p) 


0 

(AT)I u 

0 


r P ( P ) 

lu 

o 


w c (p, k) 



+ 


"7wl(p,fc)‘ 
7tu2(p, k) 
.7»3(P, k ). 


w*(p, k) 


(9a) 


( y P (p>fc) 
< yu(p,*) 
L 3/z(p,^) 



(9b) 


in which & p and T p are discrete matrices correspond- 
ing to A p and B p , respectively, j wi ( i = 1, 2, or 3) 
represents discrete matrices for the process noise, AT 
is the discrete sampling period (which is 0.0125 sec 
for this design), integer k represents the present sam- 
ple time, and the other variables are the same as 
previously defined. A perfect discrete transforma- 
tion from equation (5a) to equation (9a) will have 
slightly different coefficients in the control matrix ele- 
ments and one term of the state matrix. (See ref. 17.) 
The terms shown in equation (9a) are incorporated 
to accurately model the hold function in the digital 
computer operation. 

The control v c is related to the outputs by the 
feedback-gain matrix as 

( y P (p>.fc) I 

v c (p,fc) = -[Ky(p) K u ( p) K z (p)}< y u {p,k) \ (10) 

l 2/*(p,fc) ) 


in which K y is the proportional gain matrix, K u is 
the control filter gain, and K z is the integrator gain. 


Further partitioning of gives three terms 

Ky(p)=[K a (p) K q ( p) K nz ( p)] (11) 

for the a, q, and n z outputs. A block diagram of the 
feedback control structure described in equations (9) 
and (10) is shown in figure 5. This controller is a 
proportional integral filter (PIF) structure (refs. 17 
and 18) that is used for design and linear analysis. 
Later in this paper, an incremental structure will be 
shown for the implementation model. The general 
dynamic form of equations (9) and (10) is 

x(p, k + 1) = $(p)x(p, k) + r(p)v c (p, k) 1 


+ r w (p)w(p,fc) 

1 

> 

(12) 

y(p, k) = C(p)x(p, k) + rj( p, k) 

J 



v c {p,k) = — K(p)y(p, k) 



(13) 


Cost Function 

An optimal design process is used in the variable- 
gain output feedback algorithm. The^ algorithm re- 
quires that a discrete cost function J^(p,K(p)) be 
defined at each operating condition as 

N 

Jd( P> K (P)) = 2(N + l ) E [ xT(p ’ k + 1 )Q(P) X (P ’ k + 

fc— 0 

+ 2x t (p, k) A m (p)v c (p, k) + v£ (p, k)Jt(p)v c (p, k)] 

(14) 

in which Q(p), A m (p), and R(p) are the discrete 
weighting matrices for the states, cross product of 
states and controls, and controls, respectively, and E 
is the expectation operator. The main objective is to 
minimize a global cost J(K(p)) expressed by 

M 

J(K(p)) = E/^(p,K( P)) (/; > 0)(15) 

J= 1 

in which each of the local costs are summed and 
weighted by fj to assign relative priorities to the M 
individual operating conditions. For the design in 
this paper, all fj = 1.0. 

The design approach is to first provide optimal 
weights at each operating condition in the continuous 
domain and then to discretize these weights (refs. 19 
and 20) for use in the design algorithm (eq. (14)). 
The discrete cost function is the intermediate step to 
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minimize a continuous cost J c ( p) that is quadratic in 
outputs and controls and has the form (ref. 9) 

J c (p) = lim -I [ f E [y T (p, t)Q(p)y(p, t) 

+ vj (p, t)Rv (p)u c (p, t)] dt (16) 

in which t is time, Q is the continuous output- 
weighting matrix, and Ry is the continuous control- 
weighting matrix, which is a scalar for the HARV 
design. The outputs are related to the states by 
equations (12). 

The continuous output-weighting matrix has di- 
agonal elements as 

Q(p) — diag[Q 2/ (p) Q u ( p) Qz(p)] (17) 

for the proportional weight, control filter weight, and 
integrator weight, respectively. The proportional 
weight is further partitioned into three terms 

Q y (p) = diag[Q a (p) Q g (p) Qn 2 (p)] (18) 

for the a, q, and n z measurements. The weights for 
all design cases will be shown later in this paper. 

Stochastic- Weighting Process 

The variable-gain algorithm is stochastic and in- 
cludes adjustments for both process noise and mea- 
surement noise. Terms in the process noise cal- 
culation can include randomly distributed initial 
condition errors, control input random pseudonoise, 
random external disturbances, and integrator noise. 
The total process noise covariance matrix W is cal- 
culated as 

W = $(p)X 0 (p)$ T (p) 

+ J Q AT {e A{p)t [b(p)Vo(p)B t (p) 

+ B(p)W 0 (p)B T (p)]e / (P) t }* (19) 

and the noise terms are assumed to be white and 
have zero mean as 

£{x 0 (i,p)} = O £{x 0 (j,p)} = O 

£{xo(i,p)xoO',p)} = Xo(p)«tf 

E{v 0 (i, p)} = 0 E{v 0 (j, p)} = 0 

E{v 0 (i, p)«jf (j, p)} = Vo(P ) S ij 



E{w 0 (i, p)} = 0 E{w 0 (j, p)} = 0 

£{w 0 (i, p)w^ ( j , p)} = Wo(p)fii, 

£{v(i,p)} = 0 £{v(j,p)} = 0 

E{v(i,p)v T (j,p)} = V(p )Sij 

in which xq and Xq are the initial condition er- 
ror state vector and state covariance matrix, re- 
spectively; uq and Vo are the control input random 
pseudonoise and covariance, respectively; wq and 
Wo are the random external disturbance noise vector 
and covariance matrix, respectively; v and V are the 
measurement noise vector and covariance matrix, re- 
spectively; and $ij is the Kronecker delta. All of the 
noise terms are assumed to be uncorrelated. 

Design Conditions 

Thirty-nine design conditions (table II) are used 
for the feedback controller: 14 conditions are at 

15 000 ft, 13 conditions are at 25 000 ft, and 12 con- 
ditions are at 35 000 ft. Nineteen of the design condi- 
tions are at 1 g (Earth axis) flight. The other cases are 
at various non-lg conditions. The parameters in ta- 
ble II include the design case, altitude, Mach number, 
a, n z , and open- loop short period frequency. Because 
n z is referenced to the Z - axis of the airplane, the n z 
value for 1 g, Earth-axis flight is equal to the cosine 
of the pitch attitude, which results in a value less 
than 1.0. Most of the non-1 g trim cases are at higher 
loads; however, design cases such as 13 and 14 are 
at lower load factors. The appendix contains linear 
open-loop airplane models for all 39 design conditions 
plus the transfer functions for actuators, sensors, and 
antialiasing filters. 

Design Weights 

A design procedure for selecting the optimal cost 
weights is described in reference 7. The optimal cost 
weights for the HARV longitudinal feedback con- 
troller design are shown in table III. The first col- 
umn contains the design case number, the next five 
columns of the table contain the continuous model 
output feedback weights at each operating condition, 
and the last column contains the continuous model 
control weights. Table IV contains the stochastic 
weights for the design. Because all 39 design con- 
ditions have the same stochastic weights, weights for 
only one model are shown. The design model has 
17 states and 5 outputs; values of these stochastic 
weights are arbitrary and most of the small num- 
bers have been incorporated to avoid numerical prob- 
lems with the design algorithm. The weight on the 
first output measurement a is an order of magnitude 
larger than the other weights because a is considered 
to be a less accurate measurement. 
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Feedback-Gain Matrix 

Formulation of the feedback-gain matrix is de- 
scribed in references 8-10; only the final equation is 
included. The feedback-gain matrix, which is used 
in equation (13), has a linear functional relationship 
with p and contains both constant- and variable-gain 
parts as 

Q 

K(p) = K 0 + (24) 

i — 1 

where the variable v represents a vector of mea- 
sured variables that the designer selects for the gain- 
schedule parameters. The relationship between the 
Pi and v may be either linear or nonlinear. 

The HARV design has six gain-schedule parame- 
ters Pi(y), These parameters are functions of a, im- 
pact pressure Q c , and static pressure P s . The pi{y) 
and their limits are as follows: 


pi = 0.1a 

(1.5 < a < 65) 

P2 = 0.01 Q c 

(10 < Q c < 470) 

p 3 = O.OOlPs 

(498 < P s < 1200) 

drltf 

II 

S 

(0.008 < pa < 0.4) 

P5 = 0.1a — 3.5 

(a > 35) 

o 

II 

d? 

(a < 35) 

PG = 0.01Q C - 2.5 (Q c > 250) 

£ 

II 

o 

(Qc < 250) 


The feedback gains change continuously with the 
measured variables; the function is completely con- 
tinuous and also smooth except at two points. The 
first four parameters are used throughout the entire 
HARV flight envelope. Parameter P5 is used only 
when a is greater than 35°; parameter p§ is used 
when Q c is greater than 250 lb/ft 2 . Both p$ and pq 
have lower limits of zero and are not differentiable 
at the break points. When any value of i ^ exceeds 
the design limit, the variable is limited to the value 
shown in equations (25). The scale factors have been 
selected to keep pi(v) near unity. 

Table V contains the values of the six gain- 
schedule parameters for the 39 design conditions. 
These values are input to the variable-gain design al- 
gorithm along with the weights described previously 
to calculate the feedback gains in equation (24). Re- 
sultant values for the K* matrices (i — 0 to 6) are 
shown in table VI. Each matrix has five components 


for which the columns are consistent with the five 
output measurements. 

Linear Analysis of Feedback Controller 

This section presents results of the linear anal- 
yses for the feedback controller. Three design con- 
ditions have been selected for example cases. The 
three examples are for design cases 15, 17, and 19 
(see table II) corresponding to a low-a, high-Mach 
case (Mach = 0.7, a — 3.58°), an intermediate-o 
case (Mach = 0.33, a — 20°), and a high-a case 
(Mach = 0.26, a = 50°), respectively. All three cases 
were chosen at 25 000 ft because altitude variation 
has negligible effect. The control structure shown in 
figure 5 is used in the linear analysis. The configura- 
tion in figure 5 is used for most linear analyses but 
the pitch thrust-vectored washout filter is included 
at the output of the PIF controller for the servo- 
elastic and control sensitivity robustness analysis. 
(See discussion in “Continuous model.”) Inclusion 
of the washout filter within the controller allows the 
single control y u to be split into two separate controls 
and the control loop to be broken at the stabilator 
6 sc and thrust vector 6 VC commands. 

Feedback Gains 

Table VII contains the feedback gains for the 
three example cases. These gains are computed from 
equation (24) with the gain-schedule parameters in 
table V and the variable-gain components in table VI. 
Figure 6 is a more comprehensive picture which 
shows the gains at 25 000 ft as a function of a . The 
solid line represents the 1 g cases and the crosses 
represent gains at the other design conditions. At 
high-a conditions, the gains are essentially linear 
with a, which indicates that a is the dominant gain- 
schedule variable for this condition. Note that the 
gain plots are highly nonlinear for the a range and 
the load factor has a major effect on the gain values. 
Additional data (not included in this paper) show 
that for the HARV altitude range of 15 000-35 000 ft, 
altitude has a minor effect. 

Closed-Loop Eigenvalues 

Closed-loop eigenvalues for the three example 
cases are shown in table VIII. Row 1 in each case cor- 
responds to the normal accelerometer filter, rows 2 
and 3 correspond to a second-order pitch rate fil- 
ter, rows 4-9 correspond to actuator dynamics, and 
row 10 in design cases 15 and 17 corresponds to an 
a filter. In design case 19, the a filter state has been 
combined with another state to form a complex pair 
(rows 10 and 11). Rows 16 and 17 in each case cor- 
respond to the closed-loop phugoid. The closed-loop 
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phngoid is stable for all 39 design conditions even 
though the open-loop phugoid is unstable at a few 
design conditions. All other closed-loop eigenvalues 
correspond to combinations of short period dynam- 
ics and controller dynamics (integrator, control posi- 
tion filter, and washout filter). The other closed-loop 
eigenvalues have good damping ratios even though 
the open-loop short period is unstable for design 
case 19. 

Stability Margins 

Single loop stability analysis was done for both 
the airplane input and outputs. The input analysis 
was done by opening the control system at the single 
control loop before the split between the stabilator 
and pitch thrust-vector washout filter. The plant 
input during this analysis was augmented by the 
washout filter. Table IX shows that for the three 
example cases, the gain and phase margins are much 
better than the design guidelines of ±6 dB and ±45°, 
respectively. (See ref. 21.) These margins exclude the 
very low-frequency range of the phugoid mode. Only 
the pitch rate loop data are shown for the output 
because analysis of individual output measurement 
loops shows pitch rate is the most critical loop. 

The analysis was made with a database of 
177 cases. This data include flight conditions at three 
altitudes (15 000, 25 000, and 35000 ft), at various 
Mach numbers (from Mach 0.9 down to Mach 0.1), 
at a’s up to 65°, and other special cases that include 
heavyweight airplane, lightweight airplane, and se- 
lected conditions at the higher altitude of 40 000 ft. 
Analysis of the 177 cases shows that the gain and 
phase margins are within the guidelines for all cases 
within the HARV flight envelope (Mach < 0.7). Re- 
sults show that the nondesign cases are equally as ro- 
bust to gain and phase changes as the 39 design cases, 
which illustrates that the controller is not tuned to 
any specific design point. 

Robustness 

Structured singular value analysis, in which pa- 
rameters are modeled as uncertain gains, is useful in 
the assessment of the stability of a linear feedback 
control system when all parameters vary indepen- 
dently over given intervals. (See ref. 22.) In this 
paper, the robustness analysis is based upon a new 
algorithm for real //. (See ref. 23.) The inverse of // 
is plotted in all figures relating to robustness because 
the minimum value is indicative of the parameter per- 
centage change for which the control system remains 
stable. 

The first check is for a multiplicative error at the 
plant output with all three output loops opened. Fig- 


ure 7 contains the inverse real-// frequency plot for 
the three example cases. For the case of Mach 0.7, 
the minimum value is slightly less than 1.0, which in- 
dicates that the control system will remain stable for 
almost a 100-percent change in all plant output gains. 
The cases with an a of 20° and 50° are the most criti- 
cal with minimum values of approximately 0.6, which 
indicate that the feedback control system will re- 
main stable for approximately 60-percent changes in 
plant output gains. Although no established guide- 
lines exist for this type of analysis, these minimum 
values are considered quite good and indicate a rea- 
sonably robust control system. The complex-// anal- 
yses of 39 design conditions had similar numerical 
results with the exception of the phugoid and also 
showed that the worst cases were at an a between 
20° and 50°. 

A real-// analysis also was used to evaluate the 
sensitivity of the airplane short period stability 
derivatives and their corresponding control deriva- 
tives that affect the a and q accelerations. In the 
analysis of the four stability derivatives (fig. 8), the 
minimum values were above 1.0 for the entire fre- 
quency range. Results for the four control derivatives 
show that all minimum values were above 0.9. (See 
fig. 9.) These results are considered excellent and 
again indicate a robust feedback control system. 

Loop Transfer and Servoelastic Frequency 

Response 

A singular value loop transfer analysis was made 
to evaluate the rigid-body control system frequency 
response and crossover frequency at each design con- 
dition and also the servoelastic frequency response 
at several design conditions. The best location to 
break the plant input loop is at the single control 
point before the loop splits into two separate con- 
trols. This approach was taken for the rigid-body 
loop transfer analysis. The servoelastic analysis was 
more complex but the combination of the servo- 
elastic and rigid-body models made it easier to open 
the loops at the two control inputs. Each analysis is 
described below. 

To conduct the rigid-body singular value loop 
transfer analysis, the control loop was opened first 
at the single plant control input and then at the 
plant output (three outputs simultaneously). Be- 
cause there is only one effective control, the Bode 
response and the singular value response are identical 
at the plant input. At both input and output loca- 
tions, the crossover frequencies for all 39 design con- 
ditions ranged between 2.4 rad/sec and 10.0 rad/sec, 
with a slope of approximately 6 dB/octave. The 
highest crossover frequency is at the highest Mach 


10 



number. To obtain lower crossover frequencies at 
higher a conditions, the feedback controller filter 
bandwidth was designed to be lower by reducing the 
filter gain. Figures 10 and 11 show the crossover 
frequencies and servoelastic response for the three 
example cases. 

The servoelastic analysis was conducted as an 
extension of the rigid-body analysis; however, the 
plant input was opened at the two controls. First, a 
servoelastic modal model was placed in parallel with 
a series combination of the rigid-body and actuator 
dynamics; then, the new model was connected in 
series with the output filters. The unaugmented 
singular value servoelastic model response for a heavy 
airplane is shown in figure 10, which shows peaks of 
approximately —18 dB at 75 rad/sec and —7 dB at 
100 rad/sec. These peaks are mainly the result of 
the transfer function (not shown) from control input 
to q output (units of sec -1 ) which exhibits peaks of 
—20 dB at 75 rad/sec and —8 dB at 100 rad/sec. In 
addition, the transfer function from control input to 
n z output (units of g/ deg) exhibits peaks of -45 dB 
at 37 rad/sec, -38 dB at 50 rad/sec, and -33 dB 
at 100 rad/sec. Singular-value control loop transfer 
analyses at both the plant input and plant output 
(figs. 10 and 11) show that the structural mode at 
100 rad/sec was attenuated to approximately -20 dB 
(worst case is -19 dB at Mach 0.7). This attenuation 
is significantly greater than the guideline of —8 dB. 
(See ref. 21.) The total attenuation is mainly due to 
a combination of actuator dynamics attenuation and 
feedback controller attenuation. 

Feed- Forward Controller Design 

Feed-Forward Command Generator 
(FFCG) 

Reference 7 contains a detailed description and 
derivation of the feed-forward control equations. 
This section summarizes the controller design and 
the HARV design changes made after reference 7 
was written. The feed-forward command generator 
(FFCG) converts the pilot stick command into an 
equivalent command y cm< \ that can be interpreted by 
the feedback controller. The FFCG selects either an 
n z c or an a c command mode and makes a smooth 
transition between the two modes without any addi- 
tional work by the pilot. 

The FFCG continually generates two solutions. 
One solution is based upon the n 2?c mode with a 
stick sensitivity of 1.0 gf in. (in ref. 7 the sensitivity 


was 1.3 g/ in.) and is generated from 
3 

Vcmd = ^—(1*0 6 sp ) + ttoc (26) 

The other solution is based upon the ol c mode with 
a stick sensitivity of 10°/in. and is generated from 

3 

E Ci 

2/crn d = (10M + 20 (27) 

The C{ s are steady-state-normalized coefficients for 
a, g, and n z , respectively, 6 sp is the pilot stick posi- 
tion, and a oc is a trim that must be estimated. The 
a c mode is biased at 20° to yield a maximum com- 
mand of a — 70° with 5 in. of pitch stick. The Q’s 
are solved with the functional relationships that are 
calculated off-line by a least-squares solution. The 
methodology for derivation of the C^’s is reported in 
reference 7. 

Implementation of the FFCG is illustrated in 
figure 12 and is modified from the implementation 
shown in reference 7. This modification eliminated 
the lockout feature that was installed to minimize a 
jump between command modes during certain flight 
conditions with negative 6 $ p . The slopes of equa- 
tions (26) and (27) are compared ( 6 sp is the same 
in each equation) and the smallest absolute value 
of A y cmd is selected for the linear part of the com- 
manded output. During low-speed, high-a flight the 
bias is 20°, but during high-speed, low-a flight the 
bias a oc must be estimated. The equations and co- 
efficients for estimating a 0 c were determined from 
an off-line least-squares solution. This estimate is a 
function of impact pressure Q c and incremental load 
factor command A n ZyC ] the estimate is separated into 
a 1 g bias term as 



and an incremental bias term A% as 

Aa oc = 6.896806Q n2 + 16.2912 Ql z - 10.09235g^ 
+ 2.0293780^ (29) 

where 

Q nz = 1 °°^ 7t2 ’ c (90 <Q C < 460; A n 2 , c > 0) (30) 

Qc 
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The 1 g and incremental bias terms are then summed 
and limited to 20°. One of the inputs A n ZyC for 
this estimated bias is passed through a 1 rad/sec 
low-pass filter, which was selected as a compromise 
between time response and noise bandwidth caused 
by the pilot moving the pitch stick. By splitting 
the bias estimate into two parts, a very accurate 
estimate is obtained at 1 g, which allows the pilot 
stick to be essentially at null position independent of 
speed. Very small changes in 6 sp will balance out any 
residuals. The solution is valid up to approximately 
a = 30°. 

Command Generator Tracker (CGT) 

The objective of the CGT is to make the regu- 
lated plant outputs optimally track the linear com- 
mand model outputs. An equation is shown in ref- 
erence 7 for the calculation of the CGT feed-forward 
gain defined therein as El. The approach was to 
calculate El at each design point and then to use 
either an interpolation method or an off-line least- 
squares solution to generate an equation for calculat- 
ing El during flight. However, El was too sensitive 
to pilot stick commands at several areas of the flight 
regime, which resulted in saturation of the controls. 
In the current design, El has a constant gain of —40. 
The resulting command is then limited to maximize 
agility by allowing the controls to just reach satura- 
tion at maximum 6 sp commands. The gain of —40 
for El is representative of the calculated gains at the 
various design conditions. Adjustment of the vari- 
able command limit by trial and error resulted in the 
special conditions shown in figure 13. 

The variable limit is mainly a function of Q c 
with various modifications. The limit is maintained 
between 10 and 60 deg/sec with the upper value of 
60 deg/sec chosen to maximize agility at Mach 0.6 
and h = 25 000 ft. The lower value of 10 deg/sec is 
chosen to ensure that some agility can be attained 
at higher-a flight. Modifications to the limits are >, 
outlined below: 

1. When static pressure P s < 786 lb/ft 2 (altitude 
h > 25 000 ft) and a < 30°, the limit is reduced 
linearly with P s , which results in a smaller CGT 
command. At 35 000 ft, the reduction factor is 
approximately 0.5 because the CGT command 
was too sensitive at higher altitudes. 

2. During any sampling period (AT = 0.0125 sec) 
when the change in filtered y cmc j is less than —1 
and a is greater than 30°, the limit is increased 
linearly with a . This ensures good agility with 
a large pitch-down command (forward d sp com- 
mand) when at high a . The limit is increased by 


a factor of 4.5 when a is 60°; a typical scenario is 
for the limit to increase from 10 to 45. This limit 
increase occurs only for large and fast forward 
stick command changes (filtered y cm d rates more 
negative than —80, which corresponds to nominal 
S sp rates of —8 in/sec). 

Nonlinear Configuration 

This section briefly describes the controller con- 
figuration and puts the various parts into per- 
spective. The incremental implementation for the 
feedback controller is shown. In addition, special 
considerations for nonlinear simulation and flight are 
discussed. 

Figure 14 shows a block diagram of the HARV 
longitudinal controller containing the main compo- 
nents described in previous sections for the feedback 
and feed-forward controller designs. Pilot commands 
S S p are input to the FFCG, which generates y cm( ± 
based upon the selection of either a c or n ZjC - Out- 
put y cm d goes to both the CGT and the variable- 
gain feedback controller comprised of the PIF con- 
trol structure and gain functional, which generate 
feedback-gain matrices K(p) as a function of parame- 
ter p. The thrust- vectoring washout filter is attached 
to the output of the PIF control structure. Control 
commands 6 SC and 6 VC are input to the plant and gen- 
erate airplane responses as measured by the sensors. 
Measured y p signals are fed back to the PIF con- 
troller while other measurements are used to calcu- 
late p for both the feedback-gain functional and the 
FFCG. The p used for the FFCG is a scaled subset of 
those used for the feedback controller. An additional 
part of the controller (not discussed previously) is the 
flap scheduler, which generates commands for both 
the leading- and trailing-edge flaps as a function of 
a and air data measurements P s and Q c . The flap 
schedule is the same as used on the F-18 airplane. 

The feedback controller uses an incremental im- 
plementation (refs. 17 and 18), which is shown in fig- 
ure 15. In this implementation, K^(p) multiplies the 
incremental change in y p , K z ( p) multiplies the differ- 
ence between the sum of the measured feedbacks and 
y C mcb an d Ku( p) is incorporated into the discrete fil- 
ter loop. One advantage of the incremental approach 
is that sensor biases are subtracted out in the pro- 
portional feedback loop. In the integrator loop, the 
pilot can move the pitch stick slightly to compensate 
for biases. Position limiters are incorporated to pre- 
vent windup in the rate-to-position integrator. The 
discrete dynamics in the 6 VC actuator loop represent 
the Tustin transformation for a low-pass filter with 
a bandwidth of 1 rad/sec. The implementation from 
y u to 6 VC represents a limited washout filter. 
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Reference 7 contains a derivation 'for nonlinear 
pitch rate compensation which has two terms: one 
term for gravity compensation and the other term for 
kinematic compensation. The kinematic compensa- 
tion has been deleted in this paper because it is a 
function of sideslip angle /?, which is not a readily 
available airplane measurement. The result of this 
deletion is a slight deterioration of regulation during 
360° stability-axis rolls, mainly at a of 45° or greater. 
A second modification in the controller is a multi- 
rate approach to reduce computation time within the 
flight computer. Variables that are being computed 
every other sample period (0.025 sec as opposed to 
0.0125 sec) include a, pi( v) (i = 1 to 6), K(p), Q 
(i = 1 , 2, or 3), and some trigonometric functions for 
the nonlinear pitch rate compensation. 

Nonlinear Batch Simulation Results 

Several nonlinear batch simulations were con- 
ducted to evaluate pitch agility, a-capture capability, 
and a regulation during 360° stability-axis rolls. 

Design Guidelines 

Several preliminary open-loop (no pilot) large 
amplitude design guidelines have been established to 
assist in the design effort. Most of the guidelines that 
relate to agility are discussed in reference 11. Agility 
guidelines include minimum achievable pitch rate q 
and pitch acceleration q criteria for both pitch-up 
and pitch-down maneuvers at an altitude of 25 000 ft. 
Pitch-up agility guidelines specify the minimum peak 
q and q values for several initial trim a’s based upon 
maximum aft stick ramp command and full-throttle 
(afterburner (A/B)) command. The stick command 
is delayed for approximately 2.0 sec after throttle 
command to allow thrust to build up. The maximum 
q should be obtained within 1.0 sec from the onset 
of the pitch stick command and the maximum q 
should be achieved within 1.75 sec. The pitch-down 
agility guideline also specifies peak q and q values 
when starting from an initial trim a of 60° with a 
full-forward stick command. Values for the various 
guidelines are given in the appropriate section. 

The open- loop a-regulation tests are used to eval- 
uate how well the longitudinal controller regulates a 
while rolling about the velocity vector with a full- 
lateral stick command. The test is made at various 
initial a- trim conditions (low- a to high-a) and guide- 
lines (ref. 11) are specified at both the 90° and 360° 
wind-axis bank angles. 

Pitch-Up Agility 

The objective of this maneuver was to evaluate 
the pitch-up response (fig. 16) to a maximum aft stick 


input of 5 ill. when trimmed at an altitude of 25 000 ft 
and Mach number of 0.6. In particular, the peak 
amplitude and corresponding time of the q and q re- 
sponses were compared with design guidelines. To 
simulate the maneuver, maximum throttle was com- 
manded at time equal to 0.01 sec; 2.0 sec later, after 
thrust had built up, the pitch stick was ramped to 
maximum within 0.3 sec to simulate the approximate 
response time of the pilot. The top plot in figure 16 
shows a reaching 60° in less than 2.0 sec and then 
slowly climbing to 70° after a slight bobble. The 
sluggish response after a reached 60° was caused by 
the saturated actuator commands. Although the sta- 
bilator remains saturated, the thrust- vectored com- 
mand came out of saturation at approximately 5.6 sec 
as the a slowly converged toward 70°. The pitch 
rate q peaked at the desired guideline of approxi- 
mately 55 deg/sec, whereas the pitch acceleration q 
was greater than the guideline of 96 deg/sec 2 . The 
FFCG started in the n z mode and made a smooth 
transition to the a mode at approximately 3.0 sec; 
based on this, the pilot is unlikely to detect the tran- 
sition. The smooth time response also indicated good 
integration between the feedback and feed-forward 
controllers. 

Figures 17 and 18 illustrate q and q agility, re- 
spectively and their guidelines (ref. 11) for various a 
at 25 000 ft. The plots show that q meets or exceeds 
the guideline until initial trim a increases to approx- 
imately 45° where q drops slightly below. Similarly, 
q meets or exceeds the guideline until initial trim a 
increases to 55° where q drops slightly below. There 
is currently no guideline for a greater than 55°. 

Pitch-Down Agility 

The objective of this maneuver was to evaluate 
the pitch-down response (fig. 19) to a full- forward 
pitch stick input (—2.5 in.) starting from an a trim 
of 60° and at an altitude of 25 000 ft. The a 
decreased to 10° in approximately 2.0 sec and crossed 
0° at approximately 2.4 sec, which is slightly slower 
than the pitch-up response. Both q and q were 
significantly better than the tactical guidelines of 
—24 deg/sec and —14.3 deg/sec 2 , respectively. (See 
ref. 11.) Based on these results, the gi response is 
shown to be well damped. 

Angle-of- Attack Capture 

Figures 20 and 21 illustrate the a captures for 
altitudes of 15 000, 25 000, and 35 000 ft. Each fig- 
ure shows the stick position and the a response for 
the three altitudes. Figure 20 is for a moderate an- 
gle change; the a starts at 20° and passes maxi- 
mum lift to approximately 45°, which is based upon 
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a stick position change from 0 to 2.5 in. The ini- 
tial transient response is very similar for all three 
altitudes and reaches the desired value in just over 
2.0 sec. At t = 10 sec (nearly steady-state), the am- 
plitude difference between the upper- and lower-a 
plots is 2.8°. This small difference in amplitude is 
due to differences in the command generator feed- 
forward command gain (eqs. (26) and (27)), which is 
estimated as a function of four parameters a, Q c > Ps, 
and Qc/Ps • A maximum pilot stick position change 
of 0.28 in. would make the final steady-state ampli- 
tude the same at all altitudes. The three time re- 
sponses have very little overshoot, which indicates a 
well-damped feedback control system at all altitudes. 

Figure 21 contains similar data for a large ampli- 
tude change from 5° to nearly 60°, which is based 
upon a stick position change from 0 to 3.85 in. The 
initial transient response is again similar for all three 
a cases but reaches a maximum in approximately 
2.5 sec, which is slightly slower than the previous 
case. When the time equals 10 sec, the a amplitude 
is 56.2° at 35 000 ft, 58.7° at 25 000 ft, and 62.3° 
at 15 000 ft. Similarly to the previous case, this 
amplitude difference is due to y cm <i (not shown but 
57.0, 58.4, and 60.8, respectively), which varies with 
the estimated feed-forward gain as well as stick posi- 
tion. This small difference between a and y cm d is due 
mainly to g, which is also regulated by the feedback 
controller. Worst case pilot stick changes of approx- 
imately 0.38 in. can make the a amplitude equal to 
60° in all cases. All time responses again illustrate a 
well-damped feedback control system. 

Figure 22 is for a pitch- down a response from an 
initial trim of 60° to a capture at a = 20°; the pitch 
stick is ramped from approximately 4 in. to 0 in. in 
0.3 sec. Again, the a response was well damped and 
similar at all three altitudes. At approximately 5 sec, 
a began to slowly decrease. The decrease results 
from the increase in speed (Q c increases), which 
causes a oc (eq. (28)) to decrease below 20°, which 
in turn causes y C md ( e 0 s - (26) and (27)) to decrease 
by the same amount. The feedback controller was 
following y cm( $ very closely. 

Angle-of- Attack Regulation 

The objective of this maneuver was to evaluate 
the a regulation during full-lateral stick stability- 
axis rolls. Figure 23 illustrates four a-trim cases: 
5°, 30°, 45°, and 60°. The dashed vertical lines 
denote the time in seconds for wind-axis bank angles 
of 180° and 360°. In all cases, the a regulation is 
considerably better than the guideline of ±6° for a 
90° roll and guideline of ±10° for a 360° roll (ref. 11), 
although the two high-a cases do not roll fast enough 


for the bank angle to reach 360° within the 10 sec 
selected for this illustration. 

Piloted Simulation Results 

Pilot-in-the-loop evaluations were conducted with 
the NASA Langley Differential Maneuvering Simu- 
lator (DMS) described earlier. The evaluations in- 
cluded a series of piloted tasks designed to test the 
longitudinal and lateral-directional control system 
throughout the flight envelope. Tasks were designed 
specifically for this effort because no concise set of 
maneuvers were available to evaluate a configura- 
tion with such a large a envelope. Only longitudi- 
nal evaluation results are presented herein. How- 
ever, some of the lateral-directional maneuvers are 
described to address coupling. The piloted evalua- 
tions included a captures, n z captures, large ampli- 
tude rolls, and gross acquisition and tracking tasks. 
Results from the piloted evaluations are presented in 
terms of Cooper-Harper ratings (CHR) (ref. 24) and 
pilot comments. The CHR scale is a numerical scale 
from 1 to 10, with 1 being the best rating and 10 
the worst. (See fig. 24.) In practice, CHR’s from 1 
to 3 are referred to as level 1, ratings from 4 to 6 are 
level 2, and ratings from 7 to 9 are level 3. 

During initial piloted evaluations of the control 
law discussed in reference 7, objectionable character- 
istics were noted because of mode transitions from 
n z to a. Based on piloted simulation results, modifi- 
cations of the reference 7 configuration were required 
primarily in mode transition from n z to a command 
and back. These configuration changes are described 
herein and resulted in transitions from the n z to a 
command modes that were nearly undetectable by 
the pilots. Only results from the final version of the 
control law to be implemented and flight tested on 
the HARV are shown. Extensive evaluations were 
conducted at 15 000, 25 000, 35 000, and 40 000 ft 
but only the results from the 25 000 ft evaluations 
are shown herein for space considerations. The re- 
sults at other altitudes were similar except for the 
effect of control power deterioration near the max- 
imum a for the higher altitudes. The decrease in 
control power was due to reduced engine thrust at 
high altitudes. Altitude effects were barely notice- 
able below a — 60°. 

Three NASA test pilots participated in this study. 
Two have extensive air combat training and experi- 
ence; all have high-performance airplane experience. 
One pilot has many years of experience with simu- 
lated high-a airplanes and the other two pilots have 
one to two years of experience with simulated high-a 
airplanes. All pilots have experience with the use 
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of simulated within-visual-range air combat scenar- 
ios of high-a airplanes versus one or two conventional 
airplanes. 

The first step for each pilot in the task evalua- 
tion procedure was to repeat each maneuver to be- 
come familiar with the required piloting techniques, 
flight condition, and configuration. Once the pilots 
felt proficient, the task was repeated a few times 
to rate it. Ratings based upon a single good or 
bad attempt were prevented with this procedure. 
The drawback is that when flying the configura- 
tion in air combat, many flight conditions are en- 
countered that require proper task execution on the 
first try. Therefore, after the evaluations were com- 
pleted with all the pilots, limited one-versus-one en- 
gagements were flown against a basic F-18 airplane. 
No problems were found in adjusting between flight 
conditions; both predictability and controllability re- 
mained good. The one-versus-one results are not pre- 
sented herein. 

Note that CHR’s from the three pilots varied by 
only one number for 76 percent of the tasks and by 
no more than two numbers for any task. This rating 
spread was considered very good and implies that the 
tasks were well defined. When wide variations in the 
ratings are received, Cooper-Harper stated that the 
problem is usually traceable to poorly defined tasks. 

Nose-Up Angle-of- Attack Capture Tasks 

The a captures were conducted from two initial 
conditions: 1 g trim at Mach 0.6 and Ig trim at 

a = 20°. The captures from Mach 0.6 provided 
high q and a during pull-ups and assessed the effects 
of mode transition from dominant n z to dominant 
a command. The captures from ol — 20° produced 
lower q and a values and the system remained in the 
a-command mode for the entire maneuver. 

For each maneuver, the configuration was 
trimmed at the initial condition. When the simula- 
tion began, the pilot selected maximum A/B (max- 
imum thrust), rolled to a bank angle </> of approx- 
imately 45°, waited 2.0 sec for the engine to reach 
maximum power, and then pulled aggressively to 
capture the target a. From the time maximum A/B 
was applied to the beginning of the maneuver, Mach 
typically increased < 0.02. A nonzero </> was used 
during the maneuvers to prevent excessive energy loss 
and the possibility of entering a tail slide during the 
captures from high speeds to high o’s. For consis- 
tency, the nonzero bank angle was also used for the 
lower speed initial condition and smaller a changes. 

The Cooper-Harper task tolerance criteria re- 
quired that the capture be within ±4° for desirable 


and ±7° for adequate. No overshoots or undershoots 
which exceeded the above criteria were allowed. The 
task tolerance guidelines were intentionally restric- 
tive to make pilot gains high, which would aggrevate 
any pilot induced oscillation (PIO) tendencies that 
might exist. The CHR’s for each pilot are shown 
in figures 25 and 26 as a function of target a. The 
stringent desirable criterion was met (CHR < 4) for 
all but two cases. One pilot gave the 45° capture 
from Mach 0.6 a CHR of 5; a different pilot gave the 
60° capture from Mach 0.6 a CHR of 5. 

The a captures at 15° and 30° were conducted 
with less than maximum aft stick to meter the pitch 
rates for these relatively small a changes. With 
full-stick inputs, an a capture without an overshoot 
outside the desirable criterion was not possible. The 
pilots considered the use of less than maximum aft 
stick a natural technique that produced good pitch 
rates and capture times. 

The a captures at 45° were accomplished with 
a technique similar to the procedure used for the 
lower cn’s. However, more aft stick and higher rates 
could be applied. For the Mach 0.6 case, one of 
the pilots considered the rates required to obtain 
desirable criterion “sluggish.” When more aggressive, 
the pilot could only do the capture within ±5°, which 
resulted in a CHR of 5 ( adequate criterion met). 

At a — 60°, very little pitch control power re- 
mains and the high rates obtained from the Mach 0.6 
pull-up made the desirable criterion difficult to meet. 
Two of the pilots were able to meet the desirable cri- 
terion consistently; the other pilot was only able to 
meet the adequate criterion. 

For the FFCG described in reference 7, the a 
captures of 30° and 45° from Mach 0.6 presented 
a problem. For a capture of a = 30°, the n z to a 
transition occurred at 30°; for a capture of a - 45°, 
the transition occurred just before 45°. The FFCG 
stick gains of reference 7 changed noticeably when 
the mode transitions occurred, which resulted in a 
commanded a change. This problem was not discov- 
ered in batch simulation runs. The new implemen- 
tation of the CGT described herein works smoothly 
and causes no problems for these tasks. Transition is 
almost undetectable by the pilot. 

Nose-Down Angle-of- Attack Capture 

Tasks 

One nose-down capture maneuver was used to 
evaluate the nose-down rate, nose-down handling 
qualities, and transition from a command to n z 
command. The configuration was trimmed at 1 g and 
a = 60°. When the task began, the pilot pushed 
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over aggressively to momentarily capture a — 10° 
and then recover to 1 g level flight at 250 KCAS. 
The CHR’s primarily reflect the a capture of 10°. 
A recovery to 1# level flight was done to assure 
that no problems occurred with the a to n z mode 
transition. This maneuver was designed to evaluate 
nose-down pointing capability from a poststall initial 
condition. However, the maneuver was not intended 
to address minimum altitude loss, which is attained, 
by a pushover to maximum lift a = 37°. 

The Cooper-Harper task tolerance criteria re- 
quired that the a capture be within ±4° for desir- 
able and ±7° for adequate .' The resulting CHR’s 
are shown in figure 27. With the FFCG from refer- 
ence 7, this task received CHR’s of 7 from all pilots 
because of some problems not detected during the 
batch simulation runs. During batch runs, the stick 
motion was exact, which resulted in good captures for 
the original system. However, during piloted simu- 
lation, the initial stick position for the capture was 
not as precise. Small errors in initial stick position 
set up large amplitude oscillations induced by mode 
transitions caused by small correcting motions of the 
stick during the capture attempt. The FFCG dis- 
cussed herein solved this problem and yielded CHR’s 
of 3 or 4. The pilots still approached the task with 
caution because the available very high nose-down 
rates were difficult to arrest within the desirable crite- 
rion. The technique that was finally settled upon in- 
corporated slightly less than maximum forward stick 
to make the capture more predictable. The resulting 
workload was high but the task was repeatable with 
the intentionally reduced nose-down rate. All of the 
pilots liked the nose-down rate capability, thought 
they could meter the rate as needed, and liked the 
reassurance that a could be rapidly decreased for a 
quick recovery from poststall if necessary. 

The transition from a to n z command during the 
acceleration to 250 KCAS caused a very slight pitch 
bobble, which the pilots did not consider to be a 
problem. 

^-Heading Capture Task 

Maneuvers were also conducted to evaluate the 
n z mode of the control system. These maneuvers 
started with the configuration trimmed at Ig and 
Mach 0.6. When the simulation started, the pilot 
rolled to an appropriate bank angle (approximately 
60° for 2 g and 85° for 3 g) and pulled to capture the 
target n z . The g and Mach number were held for 90° 
of heading change and followed by a heading capture. 
The pilots were requested to rate the longitudinal and 
lateral-directional task separately. Only longitudinal 
results are given here. 


Cooper-Harper longitudinal task tolerance crite- 
ria required that the g capture be made and the g 
held within ±0.2 g for desirable and ±0.3# for ade- 
quate. As shown in figure 28, the CHR’s given for 
both the 2 g and 3# captures were 4 or better, which 
indicated desirable criterion was met in all cases. 
Two pilots had generally favorable comments about 
the g captures; one of the two pilots noted a slight 
overshoot tendency. One pilot commented that there 
was a tendency to overshoot and to drift away from 
the target g but did not consider this a serious prob- 
lem, just somewhat annoying. The overshoot and 
drift tendency may not exist in the g environment of 
the real airplane when compared with the fixed base 
simulator. One pilot had difficulty deciding between 
a 3 and 4 CHR for the 3 g task. The pilot considered 
handling qualities very good but thought the system 
could be improved. 

Coupling-During-Roll Evaluation Tasks 

Two sets of roll maneuvers were conducted to 
evaluate the lateral-directional control law (not dis- 
cussed herein) and coupling to the longitudinal axis 
during rolls. Changes in a during rolls can signif- 
icantly alter airplane roll performance and energy 
state as well as lead to departures (airplane not re- 
sponsive to pilot commands). Also, during gross ac- 
quisition tasks (pointing) at high ai, both the lateral- 
directional and longitudinal axes directly affect the 
capture. To augment target acquisition and prevent 
departures, a must be controllable. Therefore, dur- 
ing the roll tasks, the pilots were required to give 
CHR’s on a regulation. 

The Cooper-Harper a-control criteria during rolls 
required that a be maintained within ±2° for desir- 
able and within ±6° for adequate . The desirable cri- 
terion is extremely stringent, which makes the pilot 
gains high and potential problems with the controller 
more noticeable. The stringent criterion would also 
produce less favorable CHR’s. The two maneuvers 
were 1# and loaded rolls. In both rolls, the pilots 
used maximum lateral stick. The 1 g rolls started 
from a trimmed condition at the desired a and rolled 
through 360° to capture wings level ((f) = 0°) after 
the roll. The loaded rolls were more difficult. They 
required an initial nonzero (f> at the desired a and 
Mach followed by a roll to capture (j) = 90° on 
opposite side. Bank angle changes from 120° to 180° 
were noted during loaded rolls. The loaded rolls at 
a = 5° and Mach 0.6 and at a = 15° and Mach 0.4 
were at or near lg to give a comparison between the 
lg and loaded maneuver types. Results from the two 
types of rolls were very similar. The CHR’s from the 
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a regulation are shown in figure 29 for the 1 g rolls 
and in figures 30 and 31 for the loaded rolls. 

For all but one pilot at one a, desirable crite- 
rion was met (CHR < 4) even with the very strin- 
gent a regulation requirement. With the exception 
of a’s of 5° and 15°, difficulty in controlling a oc- 
curred only at the beginning of the rolls where small 
excursions in a were noted. These excursions were 
caused by slight longitudinal stick inputs associated 
with the lateral stick input and kinematic coupling 
(a changes which resulted from small (3 changes at 
the roll onset). Roll rate is automatically limited by 
the lateral-directional control law to prevent inertial 
coupling pitchout. Therefore, pitch surface satura- 
tion was not seen and regulation of a was very good 
with a low workload once the roll was established for 
a’s > 15°. 

At a = 5°, the roll occurred so fast that the pilots 
had no time to react to a changes. However, if there 
was no longitudinal stick input associated with the 
lateral input, a did not change much, which resulted 
in good CHR’s. 

The rolls conducted at a = 15° produced the 
highest workload, particularly for the 1 g rolls. At 
a = 15°, changes in Mach number during the rolls 
produced changes in a since the control law was in 
the n z mode. During the loaded rolls, these changes 
meant that small pilot stick corrections were required 
as the Mach number changed. During the 1 g rolls at 
a = 15°, small stick inputs caused mode transitions 
back and forth between the a and n z commands, 
which made the task difficult. However, most pilots 
were consistently able to maintain the very strict 
desirable a-control criterion during the rolls. 

Overall, the pilots commented that the ability to 
control a during rolls exceeded their expectations 
and was superior to anything they had previously 
flown. The task required very little attention except 
at the very beginning of the roll and at a = 15°. 

Tracking Evaluation Tasks 

Two tracking tasks were performed using a prere- 
corded target. The tasks were both done at approxi- 
mately 3#, with one task at Mach 0.45 (a’s from 15° 
to 25°) and one at Mach 0.6 (a’s from 10° to 20°). 
These tasks provided an evaluation of the acquisition 
and tracking capability of the configuration in the 
moderate-a range. During the low-speed (Mach 0.45) 
task, the effect of mode transitions between the a 
and n z commands on tracking were evaluated. The 
target was set 600 ft directly ahead of the test config- 
urations at the start of the run. When the maneuver 
started, the target rolled into a 3 g turn while holding 


a constant Mach number. Thirty seconds after the 
run began, the target reversed to a 3 g turn in the op- 
posite direction. From there, reversals were done at 
40 and 50 sec and the run ended at 70 sec. For the 
run to be considered valid, the range to the target 
could not exceed 1800 ft at any time and the Mach 
number could not vary by more than ±0.05. 

The Cooper-Harper task criteria were to keep the 
target within a 12.5-mrad-diameter pipper 50 percent 
of the time for desirable , and 10 percent of the 
time for adequate tracking. The tracking criteria 
were not to include the time during reversals. The 
longitudinal axis CHR’s for both Mach numbers are 
shown in figure 32. Desirable criterion was achieved 
by all pilots for both tasks. 

When these same tasks were performed with the 
control system described in reference 7, ratings were 
in the 5 to 6 region for the Mach 0.45 case due to 
PIO’s induced by mode transitions. With the current 
FFCG implementation, the pilots were unable to 
detect mode transitions. However, level 2 ratings 
were still given by two of the three pilots because 
of a tendency for the nose to drift above the target 
while tracking in the n z mode; the tendency for the 
nose to drift increased the pilot workload. The pilot 
that gave a level 1 rating did not consider the gradual 
correction for the nose-up drift to be a problem. The 
nose-up drift was caused by the change in a oc in 
response to a change in the n z command. Without 
a change in a oc , undetectable transitions between 
modes would not be possible. 

The Mach 0.6 tracking tasks did not include mode 
transitions and the pilot comments were good overall. 
However, the configuration was considered a little 
sensitive in pitch and the nose-up drift tendency was 
again noted. One pilot gave a CHR of 4 because 
a slight PIO tendency increased the workload, which 
caused a reaction to “consciously get out of the loop.” 

Departure Resistance 

The pilots attempted to depart the airplane with 
a series of predetermined inputs and other inputs of 
their choice. The system was proven to be extremely 
departure resistant. The only resemblance to depar- 
ture occurred during high-rate rolls at either low or 
high a with attempts to make large changes in a ei- 
ther with or without a reversal of the roll rate. For 
coordinated rolls (by neglecting gravity effects and 
setting (3 = 0), the relationship between body-axis 
roll rate p and body-axis yaw rate r is 

r—p tana (31) 
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Therefore, when rapid changes in a occurred dur- 
ing rolls, p and r changes are required to maintain 
coordination. The changes could not be done fast 
enough, which resulted in large a and (3 excursions 
due to kinematic coupling. Recovery from these ex- 
cursions was always immediate. 

Conclusions 

This paper describes the design, analysis, and 
simulation of a longitudinal controller for the F-18 
HARV. The following list contains the main conclu- 
sions and additional explanation where required: 

1. The output feedback gains are highly nonlinear 
for the range of angle of attack a. Load factor n z 
has a major effect on the gains, whereas altitude 
has a minor effect for the HARV altitude range 
of 15 000-35 000 ft. 

2. The closed- loop control system has good damping 
characteristics (as illustrated in the three exam- 
ple cases) even though the high-a example de- 
sign case has an unstable open-loop short period. 
Also, the closed-loop phugoid is stable for all 
39 design conditions even though the open-loop 
phugoid is unstable at a few design conditions. 

3. The closed- loop control system meets the gain 
and phase margin guidelines of ±6 dB and ±45°, 
respectively, at both the plant input and output 
for the entire HARV flight regime. Results from 
177 flight conditions show that nondesign cases 
are as equally robust to gain and phase changes 
as the 39 design cases, which illustrates that the 
controller is not sensitive to off-design conditions. 

4. The closed-loop control system has good robust- 
ness to changes in plant output gains and changes 
to stability and control derivatives. A real-p, anal- 
ysis for a multiplicative error at the plant output 
with all three output loops opened shows that 
gain changes of 60 percent are needed to drive 
the closed-loop system unstable. For the p anal- 
ysis of four stability derivatives that affect the 
short period, gain changes greater than 100 per- 
cent are required for instability; results for four 
control derivatives show that gain changes of at 
least 90 percent are required for instability. Al- 
though no established guidelines exist for this 
type of analysis, all of these results are consid- 
ered quite good and indicate a reasonably robust 
control system. 

5. The controller has good structural mode attenua- 
tion and greatly exceeds the guideline. Singular- 
value loop transfer analysis at both the plant in- 
put and plant output showed that a structural 


mode at 100 rad/sec was attenuated to approxi- 
mately -20 dB. This attenuation is significantly 
greater than the guideline of —8 dB. 

6. Nonlinear batch simulation analysis shows that 
the controller exhibits good agility and meets 
most of the guidelines for the entire a range. 
Pitch-up agility results show that q meets or ex- 
ceeds the guidelines until initial trim a increases 
to approximately 45° where q drops slightly be- 
low the guideline. Similarly, q meets or exceeds 
the guideline until initial trim a increases to 55° 
where q drops slightly below. There is currently 
no guideline above a = 55°. In addition, pitch- 
down results show that both q and q are signifi- 
cantly better than their tactical guidelines. 

7. During the pitch-up agility maneuver, the con- 
trol system will make a smooth transition from 
an n^-command mode to an a-command mode. 
The smooth time response indicates good inte- 
gration between the feedback and feed- forward 
controllers. 

8. Based upon nonlinear batch simulation analysis, 
the a regulation during full-lateral stick stability- 
axis rolls is considerably better than the ±6° 
guideline for a 90° roll and ±10° guideline for a 
360° roll, which indicates a tight control system. 

9. Nonlinear batch simulation analysis shows that 
the controller has good a-capture capability and 
consistency for the design altitude range. Sev- 
eral a captures illustrate similar transient re- 
sponses at all altitudes within the HARV flight 
envelope. Small differences in steady-state values 
result from variations in the feed- forward gains 
within the command generator; these differences 
can be corrected by small changes in pilot stick 
position. All plots show very little overshoot, 
which indicates a highly damped feedback con- 
trol system at all altitudes. 

10. Piloted simulation uncovered problems with an 
earlier control law design. The modified con- 
trol law discussed herein receives good Cooper- 
Harper ratings (CHR) with desirable criteria 
pilot responses for most tasks. The flying qual- 
ities are mostly level 1 or close to the level 1 to 
level 2 boundary. Piloted simulation shows the 
following: 

A. Almost all (19 of 21) CHR’s for the nose-up a 
captures meet the stringent desirable criterion 
of ±4°, which is considered good for this task. 

B. All pilots like the nose-down pitch rate ca- 
pability and the reassurance that a can be 
rapidly decreased for a quick recovery from 
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poststall if necessary. Maximum pitch rate 
makes meeting the desirable capture criterion 
of ±4° difficult; desirable criterion is attain- 
able by using less than the maximum rate. 

C. Cooper-Harper ratings for both the 2 g and 3 g 
captures show that desirable criterion was met 
by all pilots. 

D. A very stringent desirable criterion of ±2° was 
used for the a regulation during full lateral 
stick rolls at several different a’s and p’s. 
The two maneuvers were 1 g rolls and loaded 
rolls. Overall, pilots commented that the 
ability to control a during rolls exceed their 
expectations and the controller is superior to 
anything they had previously flown. This is 
reflected by the CHR’s of mostly level 1. 


E. Tracking evaluation tasks have CHR’s of ei- 
ther 3 or 4. Pilot comments for the Mach 0.6 
tracking tasks were good overall but the con- 
figuration was considered a little sensitive in 
pitch. During the Mach 0.45 task, there is a 
tendency for the nose to drift above the tar- 
get while tracking in the n z mode, which re- 
sulted in level 2 ratings because of increased 
pilot workload. The nose-up drift is caused 
by a change in a bias that is estimated in the 
feed-forward command generator. 

F. With this control system, the airplane is ex- 
tremely departure resistant. 

NASA Langley Research Center 
Hampton, VA 23681-0001 
May 25, 1994 



Appendix 
Open-Loop Models 

This appendix contains the 39 longitudinal open-loop airplane models for the operating conditions described 
in table II. In addition, transfer functions for the longitudinal actuators, sensor, and antialiasing filters are 
included. 

Open-Loop Airplane Models 

The following are state-transition matrices A and the control matrices B for the 39 airplane models. Each 
model has four states and two controls in the following order: 

1. V total velocity, ft/sec 

2. a angle of attack, deg 

3. q pitch rate, deg/sec 

4. 9 pitch attitude, deg 

5. 8 S symmetric stabilator control, deg 

6. 8 V pitch thrust vector control, deg 

All numbers smaller than 1.0E— 10 are shown as 0. Models 7, 20, and 33 have the 8 S control elements set to 
zero because 8 S is in saturation at those operating conditions. 


A1 

— 1.3820E— 02 
— 7.8960E— 03 
— 5.4849E— 03 
0 

B1 

— 2.0803E— 01 
— 1.7612E— 01 
— 1.2839E+01 
0 

A2 

— 1.4508E— 02 
— 9.9970E— 03 
— 3.1497E— 03 
0 

B2 

— 1.5206E— 01 
— 1.5669E— 01 
— 9.4529E+00 
0 

A3 

— 1.6234E— 02 
— 1.5100E— 02 
— 1.5360E— 02 
0 


8.7773E— 02 
— 1.0817E+00 
— 7.5943E+00 
0 

2.0660E— 03 
9.8815E— 01 
— 5.9355E— 01 
1.0 

— 5.6184E— 01 
0 
0 
0 

— 9.2081E— 03 
— 1.7350E— 02 
— 1.5944E+00 
0 



9.5986E— 03 
— 8.7803E— 0^. 
— 4.5764E+00 
0 

1.4572E— 03 
9.8797E— 01 
— 4.7959E— 01 
1.0000E+00 

— 5.6184E— 01 
0 
0 
0 

— 9.6873E— 04 
— 1.6548E— 02 
— 1.3096E+00 
0 



— 1.8176E— 01 
— 6.7351E— 01 
— 1.8918E+00 
0 

— 1.2164E— 05 
9.8873E— 01 
— 3.6678E— 01 
1.0000E+00 

— 5.6184E— 01 
0 
0 
0 
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— 1.3552E— 01 
— 1.2994E— 01 
— 6.1738E+00 
0 

1.0202E— 02 
— 2.2488E— 02 
— 1.4189E+00 
0 



— 5.6223E— 02 
— 4.2307E— 02 
— 4.3225E— 03 
0 

— 3.5131E— 01 
— 2.4410E— 01 
— 3.5517E— 01 
0 

6.8777E— 07 
9.9043E— 01 
— 1.8620E— 01 
1.0000E+00 

— 5.6184E— 01 
5.7095E— 09 
0 
0 

— 1.3776E— 01 
— 5.4475E— 02 
— 1.4900E+00 
0 

— 9.1970E— 02 
— 5.2898E— 02 
— 1.9423E+00 
0 



— 1.1579E— 01 
— 5.5977E— 02 
— 7.9017E— 03 
0 

— 2.2796E— 01 
— 1.0120E— 01 
— 3.3627E— 01 
0 

0 

9.8545E— 01 
— 3.1355E— 01 
1.0000E+00 

— 5.6184E— 01 
1.0403E— 08 
0 
0 

— 1.3844E— 01 
— 4.1732E— 02 
— 8.3743E— 01 
0 

— 1.2089E— 01 
— 5.5345E— 02 
— 1.6899E+00 
0 



— 1.7673E— 01 
— 4.6829E— 02 
— 1.4699E— 02 
0 

— 8.4210E— 02 
— 1.2849E— 02 
— 7.2936E— 01 
0 

0 

9.8455E— 01 
6.3096E— 02 
1.0000E+00 

— 5.4781E— 01 
3.4351E— 02 
— 1.8667E— 08 
0 

— 1.6849E— 01 
— 3.5964E— 02 
— 9.0676E— 01 
0 

— 1.6234E— 01 
— 4.4747E— 02 
— 1.5776E+00 
0 



— 2.1712E— 01 
— 3.3469E— 02 
— 1.8319E— 01 
0 

1.7024E— 02 
1.9644E— 03 
— 1.6056E+00 
0 

0 

9.9217E— 01 
— 1.4794E— 01 
1.0000E+00 

— 4.8505E— 01 
7.1098E— 02 
0 
0 



-1.5970E— 01 
-6.1952E— 02 
-1.6731E+00 
0 


— 1.6154E— 01 
— 1.8867E— 02 
— 1.2463E— 01 
0 


-5.3571E+00 
-4.8992E— 01 
-3.0092E+00 
0 


1.1041E— 05 
9.9270E— 01 
-5.6551E— 01 
1.0000E+00 


— 7.4320E— 01 
— 1.3302E— 01 
— 1.1481E+01 
0 


-8.5983E— 02 
-1.7313E— 02 
-1.6478E+00 
0 


— 1.3529E— 01 
— 2.3634E— 02 
— 5.6183E— 02 
0 


-3.7725E+00 
-4.5269E— 01 
-2.5589E+00 
0 


8.4638E— 06 
9.9120E— 01 
-5.1387E— 01 
1.0000E+00 


— 5.5771E— 01 
— 1.2017E— 01 
— 8.2145E+00 
0 


-8.2526E— 02 
-1.9626E— 02 
-1.6084E+00 
0 


— 3.5560E— 01 
— 2.6876E— 02 
— 1.1284E— 01 
0 


-5.2092E+00 
-1.2073E— 01 
-3.0182E+00 
0 


9.2728E— 06 
9.9055E— 01 
-7.2036E— 01 
1.0000E+00 


— 1.1147E+00 
— 1.2122E— 01 
— 7.289E+00 
0 


-1.3829E— 01 
-1.8472E— 02 
-1.5717E+00 
0 


-8.3637E— 02 
-2.9880E— 02 
-9.8652E— 03 
0 


-1.4968E+00 
-3.7019E— 01 
-1.0042E+00 
0 


3.0657E— 06 
9.9074E— 01 
-3.0333E— 01 
1.0000E+00 


22 



Bll 


2.8741E— 01 
8.1484E— 02 
3.4843E+00 
0 


— 8.9009E— 02 
— 3.3312E— 02 
— 1.8354E+00 
0 


A12 

— 2.3270E— 01 
— 3.2597E— 02 
— 2.1477E— 03 
0 


— 2.1267E+00 
— 1.4597E— 01 
-1.1324E+00 
0 


3.4749E— 06 
9.8796E— 01 
5.4560E— 01 
1.0000E+00 


— 4.6023E— 01 
4.3511E— 02 
— 4.2250E— 03 
0 


B12 

— 5.1037E— 01 
— 8.1249E— 02 
— 3.2163E+00 
0 


— 1.4977E— 01 
— 3.0518E— 02 
— 1.7783E+00 
0 


A13 

— 3.3181E— 03 
— 2.4662E— 02 
— 8.4316E— 03 
0 


2.8017E— 01 
— 4.4944E— 01 
— 6.3550E— 01 
0 


5.0069E— 04 
9.8834E— 01 
2.4997E— 01 
1.0000E+00 


— 5.5970E— 01 
8.7995E— 03 
— 6.5673E— 04 
0 


B13 

— 5.5619E— 02 
— 8.0847E— 02 
— 2.4145E+00 
0 


— 1.7506E— 02 
— 4.5433E— 02 
— 1.8266E+00 
0 


A14 

— 8.5408E— 02 
— 7.7199E— 02 
— 2.7813E— 02 
0 


2.0597E— 02 
— 2.8976E— 01 
— 1.6733E— 01 
0 


0 

9.7760E— 01 
3.7533E— 03 
1.0000E+00 


— 3.9728E— 01 
2.1499E— 01 
— 1.5131E— 03 
0 


B14 

— 2.8710E— 02 
— 9.9186E— 03 
— 1.1328E— 01 
0 


— 1.3626E— 01 
— 8.8371E— 02 
— 1.5340E+00 
0 


A15 

— 9.7850E— 03 
— 8.5016E— 03 
— 1.0039E— 02 
0 


— 5.7081E— 02 
— 7.3240E— 01 
— 4.5444E+00 
0 


1.5858E— 04 
9.9172E— 01 
— 4.0645E— 01 
1.0000E+00 


— 5.6184E— 01 
0 
0 
0 



B15 


— 1.4961E— 01 
— 1.2498E— 01 
— 8.7086E+00 
0 


— 1.4081E— 02 
— 1.8662E— 02 
— 1.6570E+00 
0 


A16 

— 1.2892E— 02 
— 1.1163E— 02 
— 1.2438E— 02 
0 


— 2.1418E— 01 
— 5.8328E— 01 
— 2.2707E+00 
0 


3.7871E— 04 
9.9204E— 01 
2.9286E— 01 
1.0000E+00 


— 5.6184E— 01 
0 
0 
0 


B16 

-1.3355E— 01 
— 1.0916E— 01 
— 6.0317E+00 
0 


1.3123E— 02 
— 2.3759E— 02 
— 1.7385E+00 
0 


A17 

— 4.7573E— 02 
— 3.0705E— 02 
— 5.961 80E— 03 
0 


— 3.4273E— 01 
— 2.0042E— 01 
— 4.1837E— 01 
0 


1.3207E— 06 
9.9325E— 01 
— 1.6415E— 01 
1.0000E+00 


— 5.6184E— 01 
— 2.5008E— 08 
1.1045E— 09 
0 


B17 

— 1.3713E— 01 
— 4.6247E— 02 
— 1.5134E+00 
0 


— 1.0277E— 01 
— 5.0505E— 02 
— 2.1825E+00 
0 


A18 

— 1.0177E— 01 
— 3.9980E— 02 
— 6.7700E— 03 
0 


— 2.4735E— 01 
— 8.6942E— 02 
— 3.5285E— 01 
0 


3.1883E— 07 
9.9013E— 01 
2.6703E— 01 
1.0000E+00 


— 5.5799E— 01 
1.4028E— 02 
— 4.8352E— 04 
0 


B18 

— 1.4826E— 01 
— 3.6851E— 02 
— 9.1317E— 01 
0 


— 1.3222E— 01 
— 4.6204E— 02 
— 1.7296E+00 
0 


A19 

— 1.5425E— 01 
— 3.4121E— 02 
— 6.7871E— 03 
0 


— 8.4298E— 02 
— 1.4329E— 02 
— 8.3303E— 01 
0 


2.0793E— 07 
9.8976E— 01 
5.6289E— 02 
1.0000E+00 


— 5.1488E— 01 
4.9099E— 02 
— 2.9045E— 08 
0 



B19 


1.9471E— 01 
3.2122E— 02 
1.0246E+00 
0 


— 1.7359E— 01 
— 3.8558E— 02 
— 1.6324E+00 
0 


A20 

— 1.9607E— 01 
— 2.4170E— 02 
— 1.8499E— 01 
0 


1.8325E— 02 
— 1.4411E— 03 
— 1.8167E+00 
0 


1.6409E— 06 
9.9445E— 01 
1.3313E— 01 
1.0000E+00 


— 4.1903E— 01 
7.4421E— 02 
0 
0 


B20 


0 — 1.3378E— 01 

0 — 6.0197E— 02 

0 — 1.5525E+00 

0 0 


A21 

— 1.0807E— 01 
— 1.5263E— 02 
— 8.8227E— 02 
0 


— 3.3520E+00 
— 3.5068E— 01 
— 2.1792E+00 
0 


2.4961E— 05 
9.9473E— 01 
— 3.9173E— 01 
1.0000E+00 


— 5.2796E— 01 
1.5453E— 02 
— 1.6769E— 03 
0 


B21 

— 4.8658E— 01 
— 9.3115E— 02 
— 7.7126E+00 
0 


— 1.1236E— 01 
— 2.3958E— 02 
— 2.1991E+00 
0 


A22 

— 9.0217E— 02 
— 1.8762E— 02 
— 4.2194E— 02 
0 


— 2.3089E+00 
— 3.2189E— 01 
— 1.8691E+00 
0 


1.0059E— 05 
9.9368E— 01 
— 3.5513E— 01 
1.0000E+00 


— 5.2796E— 01 
1.8015E— 02 
— 1.8242E— 03 
0 


B22 

— 3.6674E— 01 
— 8.3857E— 02 
— 5.4902E+00 
0 


— 1.0256E— 01 
— 2.7735E— 02 
— 2.1717E+00 
0 


A23 

— 2.4062E— 01 
— 2.1200E— 02 
— 8.3068E— 02 
0 


— 3.2907E+00 
— 9.6559E— 02 
— 2.0079E+00 
0 


1.0574E— 05 
9.9322E— 01 
— 4.9663E— 01 
1.0000E+00 


— 4.6023E— 01 
3.0214E— 02 
— 3.9556E— 03 
0 



B23 


7.3373E— 01 
8.4626E— 02 
4.8796E+00 
0 


— 1.7643E— 01 
— 2.5778E— 02 
— 2.1050E+00 
0 


A24 

— 5.5589E— 02 
— 2.4970E— 02 
— 8.0941E— 03 
0 


— 8.0799E— 01 
— 2.6540E— 01 
— 7.0241E— 01 
0 


3.0757E— 06 
9.9336E— 01 
2.0921E— 01 
1.0000E+00 


— 5.2796E— 01 
2.7017E— 02 
— 1.5478E— 03 
0 


B24 

— 1.8889E— 01 
— 5.6406E— 02 
— 2.3135E+00 
0 


— 9.0045E— 02 
— 3.8613E— 02 
— 2.0235E+00 
0 


A25 

— 1.5769E— 01 
— 2.6517E— 02 
— 4.9043E— 03 
0 


— 1.2464E+00 
— 1.1676E— 01 
— 7.5576E— 01 
0 


3.3212E— 06 
9.9138E— 01 
— 3.7520E— 01 
1.0000E+00 


— 4.6023E— 01 
4.5323E— 02 
— 3.0141E— 03 
0 


B25 

— 3.3591E— 01 
— 5.6241E— 02 
— 2.1366E+00 
0 


1.5702E— 01 
— 3.5324E— 02 
— 1.9487E+00 
0 


A26 

— 1.4799E— 03 
— 2.4318E— 02 
— 5.7837E— 03 
0 


3.7474E— 01 
— 3.1274E— 01 
— 4.2394E— 01 
0 


3.4523E— 04 
9.9165E— 01 
— 1.7267E— 01 
1.0000E+00 


— 5.5970E— 01 
9.1706E— 03 
— 4.7388E— 04 
0 


B26 

— 3.7089E— 02 
— 5.5740E— 02 
-1.5990E+00 
0 


— 1.2182E— 02 
— 4.8642E— 02 
-1.8719E+00 
0 


A27 

— 5.8649E— 02 
— 9.7055E— 02 
— 2.0472E— 02 
0 


1.4562E— 01 
— 2.7542E— 01 
— 1.0391E— 01 
0 


0 

9.8401E— 01 
2.3766E— 03 
1.0000E+00 


— 3.9728E— 01 
2.2381E— 01 
— 1.0788E— 03 
0 



B27 


1.8896E— 02 
6.8075E— 03 
7.4668E— 02 
0 


— 1.3304E— 01 
— 9.1693E— 02 
— 1.5150E+00 
0 


A28 

— 8.8257E— 03 
— 9.0749E— 03 
— 1.8778E— 02 
0 


— 2.7677E— 01 
— 4.7549E— 01 
— 2.3763E+00 
0 


— 9.4246E— 04 
9.9486E— 01 
— 2.4379E— 01 
1.0000E+00 


— 5.6184E— 01 
0 
0 
0 


B28 

— 1.3431E— 01 
— 8.3626E— 02 
— 5.4646E+00 
0 


9.7489E— 03 
— 2.1907E— 02 
— 1.8154E+00 
0 


A29 

— 1.3503E— 02 
— 1.1684E— 02 
— 1.5311E— 02 
0 


3.4369E— 01 
3.8144E— 01 
1.1206E+00 
0 


— 1.3901E— 03 
9.9544E— 01 
— 1.4936E— 01 
1.0000E+00 


— 5.6184E— 01 
0 
0 
0 


B29 

— 1.3362E— 01 
— 7.4738E— 02 
— 3.8958E+00 
0 


2.3442E— 03 
— 2.7580E— 02 
— 1.9604E+00 
0 


A30 

— 3.6378E— 02 
— 2.1764E— 02 
— 9.9397E— 03 
0 


3.3058E— 01 
1.6148E— 01 
5.0165E— 01 
0 


1.3543E— 06 
9.9540E— 01 
— 1.4496E— 01 
1.0000E+00 


— 5.6178E— 01 
1.1630E— 03 
— 4.5991E— 05 
0 


B30 

— 1.2523E— 01 
— 3.8871E— 02 
— 1.5783E+00 
0 


— 1.1121E— 01 
— 4.1541E— 02 
— 2.1698E+00 
0 


A31 

— 8.8703E— 02 -2.6830E-01 7.5841E-07 

— 2.7580E— 02 -7.2754E-02 9.9365E-01 

— 3.3416E— 03 — 3.6790E— 01 -2.2129E-01 

0 0 1.0000E+00 


— 5.4041E— 01 
2.6327E— 02 
— 9.0744E— 04 
0 



B31 


1.5815E— 01 
3.2152E— 02 
9.9877E— 01 
0 


— 1.7097E— 01 
— 4.3409E— 02 
— 2.0099E+00 
0 


A32 

— 1.3235E— 01 
— 2.3482E— 02 
— 3.4882E— 03 
0 


— 8.0098E— 02 
— 1 3921E— 02 
— 9.0622E— 01 
0 


7.6752E— 07 
9.9361E— 01 
4.7900E— 02 
1.0000E+00 


— 4.7441E— 01 
5.2846E— 02 
— 3.3557E— 08 
0 


B32 

— 2.1363E— 01 
— 2.7498E— 02 
-1.0958E+00 
0 


— 2.1852E— 01 
— 3.6346E— 02 
— 1.8787E+00 
0 


A33 

— 1.5347E— 01 
— 1.8737E— 02 
— 7.6505E— 02 
0 


1.2320E— 02 
9.2665E— 04 
— 2.4306E+00 
0 


— 7.5442E— 07 
9.9604E— 01 
— 4.2878E— 02 
1.0000E+00 


— 4.0696E— 01 
6.4489E— 02 
0 
0 


B33 


0 — 1.0623E— 01 

0 — 3.6799E— 02 

0 — 1.2328E+00 

0 0 


A34 

— 6.7577E— 02 
— 1.2335E— 02 
— 5.8820E— 02 
0 


— 1.9369E+00 
— 2.4196E— 01 
— 1.4745E+00 
0 


1.7594E— 05 
9.9631E— 01 
— 2.6241E— 01 
1.0000E+00 


— 5.2796E— 01 
1.6141E— 02 
— 1.1617E— 03 
0 


B34 

— 3.0812E— 01 
— 6.2776E— 02 
— 4.9817E+00 
0 


— 1.1210E— 01 
— 2.5575E— 02 
— 2.2365E+00 
0 


A35 

— 5.7268E— 02 
— 1.5078E— 02 
— 2.8420E— 02 
0 


— 1.2721E+00 
— 2.2124E— 01 
— 1.2731E+00 
0 


6.2707E— 06 
9.9559E— 01 
— 2.3747E— 01 
1.0000E+00 


— 5.2796E— 01 
1.8821E— 02 
— 1.2640E— 03 
0 



B35 


2.3273E— 01 
5.6403E— 02 
3.5338E+00 
0 


— 1.1296E— 01 
— 2.9213E— 02 
— 2.2052E+00 
0 


A36 

— 1.5707E— 01 
— 1.6732E— 02 
— 5.7006E— 02 
0 


— 1.9160E+00 
— 7.5129E— 02 
— 1.3094E+00 
0 


6.5864E— 06 
9.9527E— 01 
3.3135E— 01 
1.0000E+00 


— 4.6023E— 01 
3.1566E— 02 
— 2.7408E— 03 
0 


B36 

— 4.6614E— 01 
— 5.6943E— 02 
— 3.145 1E+00 
0 


— 1.8737E— 01 
— 2.7026E— 02 
— 2.1272E+00 
0 


A37 

— 7.7855E— 04 
— 1.5770E— 02 
— 4.6668E— 03 
0 


3.4732E— 01 
— 2.6909E— 01 
— 4.4189E— 01 
0 


1.3703E— 04 
9.9428E— 01 
1.4684E— 01 
1.0000E+00 


— 5.5970E— 01 
7.1889E— 03 
— 3.0138E— 04 
0 


B37 

— 4.0356E— 02 
— 4.9997E— 02 
— 1.8130E+00 
0 


— 3.1748E— 02 
— 4.4093E— 02 
— 2.1694E+00 
0 


A38 

— 1.5736E— 01 
— 1.7558E— 02 
— 3.0896E— 03 
0 


— 3.7928E— 01 
— 1.2968E— 02 
— 1.3071E+00 
0 


1.3581E— 06 
9.9421E— 01 
5.7023E— 02 
1.0000E+00 


— 3.6114E— 01 
6.3349E— 02 
0 
0 


B38 

— 3.0976E— 01 
— 3.2469E— 02 
— 1.5525E+00 
0 


— 2.3194E— 01 
— 3.2199E— 02 
— 1.9473E+00 
0 


A39 

— 7.1927E— 02 
— 3.9744E— 02 
— 5.2789E— 03 
0 


8.7245E— 02 
— 1.1375E— 01 
— 2.9716E— 01 
0 


0 

9.8892E— 01 
2.9038E— 03 
1.0000E+00 


— 3.9728E— 01 
1.1686E— 01 
— 7.4713E— 04 
0 


B39 

— 6.3100E— 02 
— 2.0587E— 02 
— 3.8539E— 01 
0 


— 1.6850E— 01 
— 5.5758E— 02 
— 1.7212E+00 
0 



Actuator, Sensor, and Filter Transfer Functions 

The stabilator transfer function is 


(36.4) 2 

(105) 2 

s 2 + 2(0.41) (36.4)s + (36.4) 2 

s 2 + 2(0.59)(105)s + (105) 2 


The pitch thrust- vectoring transfer function is 


(75) 2 

s 2 + 2(0.59)(75)s + (75) 2 


The angle-of-attack probe transfer function is 

14 

7+14 


The normal acceleration n z state-space output equation is 


n z = 


V. 

-(q - a) + 



7 r 

180 


where true airspeed V can be calculated from the Mach number (table II) and the speed of sound at the 
appropriate altitude, g is gravity (also a function of altitude), and X acc is the distance of the sensor (12.35 ft) 
forward of the center of gravity. The d and q equations are calculated from the second and third rows of the 
appropriate and matrices. 


Antialiasing filter transfer functions include the following: 


1. Angle of attack: 


(209) 2 

s 2 +2(0.74) (209)s+(209) 2 


2. Pitch rate: 


(78. 5) 2 

s 2 +2(0.89)(78.5)s+(78.5) 2 


3. Load factor: 


(200) 2 

s 2 +2(0.89)(200)s+(200)^ 


For control system analysis, all the transfer functions were included; for design, the a antialiasing filter was 
not included and the load factor antialiasing filter was approximated by a first-order filter with a 200-rad/sec 
bandwidth. The latter generates an extra state to eliminate the D i terms in the n z output equation because 
the design algorithm cannot accept terms. 
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Table I. Mass Properties 


State 

Weight, 

lbf 

x cg, 

percent c 

z cg, 

relative to WL 

Ixx-) 

slug-ft 2 

hy 

slug-ft 2 

Izzi 

slug-ft 2 

Ixz j 

slug-ft 2 

Light 

31618 

26.6 

103.4 

22 163 

172 238 

186823 

-2043 

Nominal 

35 765 

23.3 

105.4 

22 633 

174246 

189336 

-2132 

Heavy 

37619 

23.3 

105.9 

22 938 

179 130 

194003 

-2507 


Table II. Design Conditions 






Short-period 

Design case 

Mach 

a, deg 

n z , g units 

frequency, rad/sec 


Altitude = 15000 ft 


1 

0.70 

2.52 

1.0 

2.7 

2 

.60 

3.37 

1.0 

2.1 

3 

.49 

5 

1.0 

1.4 

4 

.27 

20 

.94 

.57 

5 

.21 

35 

.82 

.54 

6 

.20 

50 

.80 

.83 

7 

.22 

65 

.82 

1.3 

8 

.70 

20 

6.3 

1.7 

9 

.60 

20 

4.9 

1.6 

10 

.60 

35 

6.9 

1.7 

11 

.40 

20 

2.1 

.97 

12 

.40 

35 

3.1 

.97 

13 

.30 

5 

.37 

.79 

14 

.10 

45 

.22 

.35 


Altitude = 25 000 ft 


15 

0.70 

3.58 

1.0 

2.1 

16 

.59 

5 

1.0 

1.5 

17 

.33 

20 

.94 

.63 

18 

.26 

35 

.88 

.56 

19 

.26 

50 

.90 

.90 

20 

.28 

65 

.92 

1.3 

21 

.70 

20 

4.2 

1.5 

22 

.60 

20 

3.2 

1.4 

23 

.60 

35 

4.5 

1.4 

24 

.40 

20 

1.4 

.82 

25 

.40 

35 

2.0 

.83 

26 

.30 

5 

.24 

.64 

27 

.10 

45 

.14 

.27 


Altitude = 35 000 ft 


28 

0.70 

5.34 

1.0 

1.5 

29 

.60 

7.24 

.99 

1.0 

30 

.41 

20 

.94 

.70 

31 

.34 

35 

.94 

.60 

32 

•34 

50 

.95 

.94 

33 

.35 

60 

.95 

1.6 

34 

.70 

20 

2.7 

1.2 

35 

.60 

20 

2.0 

1.2 

36 

.60 

35 

2.9 

1.1 

37 

.40 

5 

.28 

.66 

38 

.40 

50 

1.4 

1.1 

39 

.20 

45 

.34 

.53 




Table III. Optimal Cost Weights 


Design case 

Qol 

Qq 

Qn,z 

Qu 

Qz 

Rv 

1 

94.6876 

91.4846 

14.7383 

2.4387 

11.8359 

0.0038 

2 

116.0329 

75.9189 

8.9587 

4.0312 

14.5041 

.0063 

3 

138.2935 

58.1897 

4.4272 

9.0703 

17.2867 

.0142 

4 

722.7099 

79.3696 

1.5625 

10.0000 

22.5847 

.0214 

5 

762.9073 

40.3961 

.3387 

10.0000 

23.8409 

.0238 

6 

800.0000 

3.6249 

.0171 

3.2258 

25.0000 

.0403 

7 

200.0657 

.8319 

.0129 

2.0000 

25.0082 

.0164 

8 

126.5197 

66.6656 

7.7252 

6.1515 

15.8150 

.0096 

9 

150.7106 

44.8945 

5.3054 

6.9444 

18.8388 

.0109 

10 

177.6644 

20.3447 

2.9014 

1.5379 

22.2081 

.0096 

11 

684.2341 

112.5526 

6.8553 

4.7236 

21.3823 

.0295 

12 

782.2292 

17.6359 

3.7769 

4.7236 

24.4447 

.0295 

13 

167.6470 

32.3296 

.9339 

1.7803 

20.9559 

.0445 

14 

803.2323 

.3282 

.0815 

18.1406 

25.1010 

.1134 

15 

132.2815 

59.7755 

8.8535 

4.0312 

16.5352 

0.0063 

16 

150.2881 

45.8563 

4.7661 

7.9012 

18.7860 

.0123 

17 

756.2476 

46.1103 

1.2841 

10.0000 

23.6327 

.0350 

18 

787.4826 

15.8949 

.2646 

10.0000 

24.6088 

.0443 

19 

799.4125 

4.2137 

.0158 

10.0000 

24.9816 

.0343 

20 

797.1861 

6.3590 

.0969 

10.0000 

12.4560 

.0164 

21 

158.9245 

37.8135 

4.1725 

7.9012 

19.8656 

.0123 

22 

174.8601 

23.4588 

2.5917 

9.0703 

21.8575 

.0142 

23 

189.7631 

9.8827 

1.2646 

2.2676 

23.7204 

.0142 

24 

743.9655 

56.5155 

3.1610 

6.6098 

23.2489 

.0413 

25 

793.8685 

8.2248 

1.5487 

6.4515 

24.8084 

.0403 

26 

183.1379 

17.3183 

.4543 

2.7127 

22.8922 

.0678 

27 

799.6596 

3.9695 

.0128 

20.0000 

24.9894 

.0953 

28 

163.5634 

32.9563 

4.3908 

31.6049 

20.4454 

0.0123 

29 

173.6166 

24.9270 

2.3669 

40.0000 

21.7021 

.0278 

30 

775.2657 

27.3099 

1.0665 

36.2812 

24.2271 

.0567 

31 

799.5079 

3.9247 

.2094 

40.0000 

24.9846 

.0772 

32 

799.9307 

3.6961 

.0152 

40.0000 

12.4989 

.0314 

33 

795.0167 

8.4284 

.1969 

40.0000 

6.2111 

.0109 

34 

180.8571 

18.1610 

1.8924 

40.0000 

22.6071 

.0193 

35 

189.2436 

10.5712 

1.0957 

40.0000 

23.6555 

.0193 

36 

196.5172 

3.9141 

.4792 

14.6924 

24.5647 

.0230 

37 

187.4423 

12.9098 

.5585 

10.2030 

23.4303 

.0638 

38 

798.2298 

5.4121 

.0001 

40.0000 

24.9447 

.0230 

39 

802.7795 

.8515 

.0110 

40.0000 

25.0869 

.0494 


Table IV. Stochastic Weights 61 


State 

xq weight 


State 

xq weight 


Output 

v weight 

1 

0.01 


10 

0 


1 

0.10 

2 

.01 


11 

.05 


2 

.01 

3 

.01 


12 

.05 


3 

.01 

4 

.01 


13 

.05 


4 

.01 

5 

0 


14 

.01 


5 

.01 

6 

0 


15 

0 




7 

0 


16 

.01 




8 

0 


17 

0 




9 

0 








“Stochastic weights for v$ and wq are set equal to 0. All 39 design cases have the same weights. 
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Table V. Gain-Schedule Design Parameters 



Design case 

Pi 

P 2 

P 3 

PA 

Pb 

P 6 

i 

0.2515 

4.6232 

1.1943 

0.3871 

0.0000 

2.1232 

2 

.3370 

3.2904 

1.1943 

.2755 

.0000 

.7904 

3 

.5000 

2.0866 

1.1943 

.1747 

.0000 

.0000 

4 

2.0000 

.6131 

1.1943 

.0513 

.0000 

.0000 

5 

3.5000 

.3581 

1.1943 

.0300 

.0000 

.0000 

6 

5.0000 

.3271 

1.1943 

.0274 

1.5000 

.0000 

7 

6.5000 

.3951 

1.1943 

.0331 

3.0000 

.0000 

8 

2.0000 

4.6232 

1.1943 

.3871 

.0000 

2.1232 

9 

2.0000 

3.2904 

1.1943 

.2755 

.0000 

.7904 

10 

3.5000 

3.2904 

1.1943 

.2755 

.0000 

.7904 

11 

2.0000 

1.3920 

1.1943 

.1166 

.0000 

.0000 

12 

3.5000 

1.3920 

1.1943 

.1166 

.0000 

.0000 

13 

.5000 

.7695 

1.1943 

.0644 

.0000 

.0000 

14 

4.5000 

.0838 

1.1943 

.0070 

1.0000 

.0000 

15 

0.3580 

3.0402 

0.7854 

0.3871 

0.0000 

0.5402 

16 

.5000 

2.0819 

.7854 

.2651 

.0000 

.0000 

17 

2.0000 

.6120 

.7854 

.0779 

.0000 

.0000 

18 

3.5000 

.3884 

.7854 

.0495 

.0000 

.0000 

19 

5.0000 

.3729 

.7854 

.0475 

1.5000 

.0000 

20 

6.5000 

.4512 

.7854 

.0575 

3.0000 

.0000 

21 

2.0000 

3.0402 

.7854 

.3871 

.0000 

.5402 

22 

2.0000 

2.1637 

.7854 

.2755 

.0000 

.0000 

23 

3.5000 

2.1637 

.7854 

.2755 

.0000 

.0000 

24 

2.0000 

.9154 

.7854 

.1166 

.0000 

.0000 

25 

3.5000 

.9154 

.7854 

.1166 

.0000 

.0000 

26 

.5000 

.5060 

.7854 

.0644 

.0000 

.0000 

27 

4.5000 

.0551 

.7854 

.0070 

1.0000 

.0000 

28 

0.5344 

1.9278 

0.4980 

0.3871 

0.0000 

0.0000 

29 

.7246 

1.3720 

.4980 

.2755 

.0000 

.0000 

30 

2.0000 

.6171 

.4980 

.1239 

.0000 

.0000 

31 

3.5000 

.4226 

.4980 

.0849 

.0000 

.0000 

32 

5.0000 

.4032 

.4980 

.0810 

1.5000 

.0000 

33 

6.0000 

.4498 

.4980 

.0903 

2.5000 

.0000 

34 

2.0000 

1.9278 

.4980 

.3871 

.0000 

.0000 

35 

2.0000 

1.3720 

.4980 

.2755 

.0000 

.0000 

36 

3.5000 

1.3720 

.4980 

.2755 

.0000 

.0000 

37 

.5000 

.5804 

.4980 

.1166 

.0000 

.0000 

38 

5.0000 

.5804 

.4980 

.1166 

1.5000 

.0000 

39 

4.5000 

.1408 

.4980 

.0283 

1.0000 

.0000 



Table VI. Variable-Gain Components 



Va 

Vq 

Un,z 

Uu 

Vz 

K 0 

-10.6285 

-25.4721 

-5.3189 

21.9340 

-30.8027 

Ki 

-1.2185 

-1.0865 

-10.2974 

-.0423 

-4.5770 

k 2 

-3.5173 

-12.7954 

-18.7768 

12.7372 

14.6450 

k 3 

4.4277 

.8839 

.2946 

-6.7833 

-16.0592 

k 4 

33.5886 

51.6908 

38.9416 

9.9597 

-37.4526 

k 5 

-1.5707 

-.6439 

50.1606 

1.7770 

19.0618 

k 6 

4.5668 

9.1496 

24.2615 

1.5378 

-8.7858 


Table VII. Feedback Gains 


Design case 

K a 

K q 

K n ,z 

K u 

K z 

15 

-2.8113 

-39.1150 

-37.6789 

60.0010 

-19.7742 

17 

-9.1233 

-30.7536 

-34.1390 

25.0931 

-46.5250 

19 

-15.3165 

-33.4935 

13.5122 

24.2836 

-34.0246 


35 





Table VIII. Closed-Loop Eigenvalues 


Imaginary 


Damping ratio 


1 

-195.5 

0 

1.0000 

2 

-106.9 

30.47 

.9617 

3 

-106.9 

-30.47 

.9619 

4 

-61.28 

86.62 

.5776 

5 

-61.28 

-86.62 

.5776 

6 

-44.47 

62.50 

.5797 

7 

-44.47 

-62.50 

.5797 

8 

-24.18 

33.38 

.5866 

9 

-24.18 

-33.38 

.5866 

10 

-13.89 

0 

1.0000 

11 

-11.49 

11.52 

.7061 

12 

-11.49 

-11.52 

.7061 

13 

-1.266 

.3584 

.9622 

14 

-1.266 

-.3584 

.9622 

15 

-.5100 

0 

1.0000 

16 

-.003157 

.04924 

.0640 

17 

-.003157 

-.04924 

.0640 

Design case 17 

i 

-199.3 

0 

1.0000 

2 

-76.52 

38.21 

.8947 

3 

-76.52 

-38.21 

.8947 

4 

-61.96 

85.09 

.5887 

5 

-61.96 

-85.09 

.5887 

6 

-43.12 

62.09 

.5705 

7 

-43.12 

-62.09 

.5705 

8 

-16.09 

33.73 

.4306 

9 

-16.09 

-33.73 

.4306 

10 

-13.27 

0 

1.0000 

11 

-8.941 

2.004 

.9758 

12 

-8.941 

-2.004 

.9758 

13 

-1.096 

1.213 

.6706 

14 

-1.096 

-1.213 

.6706 

15 

-.3913 

0 

1.0000 

16 

-.01082 

.1175 

.0917 

17 

-.01082 

-.1175 

.0917 

Design case 19 

i 

-200.1 

0 

1.0000 

2 

-75.29 

37.82 

.8936 

3 

-75.29 

, -37.82 

.8936 

4 

-62.00 

85.07 

.5890 

5 

-62.00 

-85.07 

.5890 

6 

-43.35 

62.10 

.5724 

7 

-43.35 

-62.10 

.5724 

8 

-15.80 

33.59 

.4258 

9 

-15.80 

-33.59 

.4258 

10 

-13.81 

1.890 

.9908 

11 

-13.81 

-1.890 

.9908 

12 

-4.417 

0 

1.0000 

13 

-.8716 

.9895 

.6610 

14 

-.8716 

-.9895 

.6610 

15 

-.3045 

0 

1.0000 

16 

-.09067 

.1351 

.5571 

17 

-.09067 

-.1351 

.5571 
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Figure 7. Results of real-fi analysis at plant output, h = 25 000 ft; 1 g trim. 
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Figure 14. Block diagram of HARV longitudinal controller. 
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Figure 21. Alpha captures of 5° to 60° for three altitudes 



Pitch stick, in. 

a, deg at 15 000 ft 

a, deg at 25 000 ft 

a, deg at 35 000 ft 








Adequacy for selected task 
or required operation* 


Yes 


Is it 

satisfactory 
without improvement? 


No 


Deficiencies 

warrant 

improvement 


Yes 


r Is adequate ' 



performance 

No 

Deficiencies 

attainable with 



require 

a tolerable 


improvement 

pilot workload? 



^Yes 


Is it 

No 

Improvement 

controllable? 

V J 


mandatory 



Aircraft 

Demands on the pilot in 

Pilot 

characteristics 

selected task or required operation* 

rating 


Excellent 
Highly desirable 

Pilot compensation not a factor 
for desired performance 

1 

Good 

Negligible deficiencies 

Pilot compensation not a factor 
for desired performance 

2 

Fair - Some mildly 
unpleasant deficiencies 

Minimal pilot compensation 
required for desired perfomance 

3 

Minor but annoying 
deficiencies 

Desired performance requires 
moderate pilot compensation 

4 

Moderately objectionable 
deficiencies 

Adequate performance requires 
considerable pilot compensation 

5 

Very objectionable but 
tolerable deficiencies 

Adequate performance requires 
extensive pilot compensation 

6 

Major deficiencies 

Adequate performance not attainable 
with maximum tolerable pilot 
compensation controllability not in 
question 

7 

Major deficiencies 

1 

Considerable pilot compensation is 
required for control 

8 

Major deficiencies 

Intense pilot compensation is 
required to retain control 

9 

Major deficiencies 

Control will be lost during some 
portion of required operation 

10 


‘Definition of required operation involves designation of flight phase 
and/or subphases with accompanying conditions 


Figure 24. Cooper-Harper handling qualities rating scale (from ref. 24). 
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Figure 25. Nose-up angle-of-attack captures from Mach 0.6 at 25 000 ft. 
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Figure 26. Nose-up angle-of-attack captures from a = 20° at 25 000 ft. 
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Figure 29. Angle-of-attack regulation during 1 g rolls at 25 000 ft. 
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Figure 30. Angle-of-attack regulation during loaded rolls at Mach 0.6 and 25 000 ft. 
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Figure 31. Angle-of-attack regulation during loaded rolls at Mach 0.4 and 25000 ft. 



Maneuver Mach number 


Figure 32. Longitudinal tracking at 3 g and 25 000 ft. 
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